cap log close
set more off

* stata-mp4-5g -b do bk-2000-resid1.do;
* stata-mp8-5g -b do bk-2000-resid1.do;

log using bk-2000-resid1.log, replace

*** Prep for Simulations
* Exported to excel to make code used below
use ethn metaread puma perwt if metaread!=0 using ./data-raw/1a-US-Master2-Wage.dta, clear
egen temp1=max(puma), by(metaread)
egen temp2=min(puma), by(metaread)
drop if temp1==temp2
ren perwt ct
save ./temp/Residential-Baseline.dta, replace
collapse (sum) ct, by(puma metaread) fast
egen temp1=total(ct), by(metaread)
replace ct=ct/temp1
drop temp1
sort metaread puma
gen sh=0
replace sh=sh[_n-1]+ct[_n-1] if metaread==metaread[_n-1]
ren puma zpuma
save ./temp/Residential-Distribution, replace
list, clean noobs

*** Real values
use ./temp/Residential-Baseline.dta, clear
collapse (sum) ct, by(ethn puma metaread) fast
egen cell=group(puma metaread)
egen msatot=sum(ct), by(metaread)
egen celltot=sum(ct), by(cell)
egen Emsatot=sum(ct), by(ethn metaread)
egen Ecelltot=sum(ct), by(ethn cell)
gen Esh=Ecelltot/Emsatot
gen Tsh=celltot/msatot
gen Edev=(Esh-Tsh)^2
collapse (sum) Edev ct, by(ethn metaread) fast
collapse (mean) Edev [fw=ct], by(ethn) fast
gsort -Edev
list in 1/10, clean noobs 
sort ethn
gen n=0
save ./temp/residential-isolation-sim1, replace

*** Randomizations
forvalues i = 1/1000 {
di `i'
set seed `i'
use ./temp/Residential-Baseline.dta, clear
drop puma
gen zpuma=.
g zsh = uniform()
replace zpuma=1000 if metaread==80 & zsh>0
replace zpuma=1101 if metaread==80 & zsh>0.2336009
replace zpuma=1102 if metaread==80 & zsh>0.380811
replace zpuma=1103 if metaread==80 & zsh>0.5528815
replace zpuma=1201 if metaread==80 & zsh>0.7011385
replace zpuma=1202 if metaread==80 & zsh>0.8561372
replace zpuma=2100 if metaread==160 & zsh>0
replace zpuma=2201 if metaread==160 & zsh>0.1762791
replace zpuma=2202 if metaread==160 & zsh>0.3078116
replace zpuma=2300 if metaread==160 & zsh>0.44504
replace zpuma=2401 if metaread==160 & zsh>0.6363766
replace zpuma=2402 if metaread==160 & zsh>0.7689905
replace zpuma=500 if metaread==200 & zsh>0
replace zpuma=601 if metaread==200 & zsh>0.2002669
replace zpuma=602 if metaread==200 & zsh>0.4200287
replace zpuma=603 if metaread==200 & zsh>0.5601357
replace zpuma=604 if metaread==200 & zsh>0.7147614
replace zpuma=605 if metaread==200 & zsh>0.8675508
replace zpuma=3600 if metaread==240 & zsh>0
replace zpuma=3701 if metaread==240 & zsh>0.1657861
replace zpuma=3702 if metaread==240 & zsh>0.4060837
replace zpuma=3801 if metaread==240 & zsh>0.576673
replace zpuma=3802 if metaread==240 & zsh>0.7701318
replace zpuma=200 if metaread==320 & zsh>0
replace zpuma=300 if metaread==320 & zsh>0.5146222
replace zpuma=101 if metaread==380 & zsh>0
replace zpuma=102 if metaread==380 & zsh>0.4183058
replace zpuma=3200 if metaread==440 & zsh>0
replace zpuma=3300 if metaread==440 & zsh>0.2607967
replace zpuma=3400 if metaread==440 & zsh>0.6893733
replace zpuma=201 if metaread==480 & zsh>0
replace zpuma=202 if metaread==480 & zsh>0.5290791
replace zpuma=700 if metaread==520 & zsh>0
replace zpuma=900 if metaread==520 & zsh>0.0384031
replace zpuma=1101 if metaread==520 & zsh>0.0636757
replace zpuma=1102 if metaread==520 & zsh>0.091776
replace zpuma=1103 if metaread==520 & zsh>0.1420465
replace zpuma=1104 if metaread==520 & zsh>0.1734739
replace zpuma=1105 if metaread==520 & zsh>0.2012936
replace zpuma=1106 if metaread==520 & zsh>0.2244878
replace zpuma=1107 if metaread==520 & zsh>0.2455292
replace zpuma=1201 if metaread==520 & zsh>0.2691039
replace zpuma=1202 if metaread==520 & zsh>0.2973433
replace zpuma=1203 if metaread==520 & zsh>0.3338203
replace zpuma=1204 if metaread==520 & zsh>0.3610452
replace zpuma=1205 if metaread==520 & zsh>0.3865996
replace zpuma=1206 if metaread==520 & zsh>0.4106683
replace zpuma=1301 if metaread==520 & zsh>0.4375068
replace zpuma=1302 if metaread==520 & zsh>0.4673585
replace zpuma=1303 if metaread==520 & zsh>0.4964662
replace zpuma=1304 if metaread==520 & zsh>0.5274713
replace zpuma=1305 if metaread==520 & zsh>0.5543562
replace zpuma=1401 if metaread==520 & zsh>0.5971929
replace zpuma=1402 if metaread==520 & zsh>0.6266309
replace zpuma=1501 if metaread==520 & zsh>0.6535808
replace zpuma=1502 if metaread==520 & zsh>0.6815675
replace zpuma=1503 if metaread==520 & zsh>0.714776
replace zpuma=1504 if metaread==520 & zsh>0.7408661
replace zpuma=1505 if metaread==520 & zsh>0.7737637
replace zpuma=1600 if metaread==520 & zsh>0.8063925
replace zpuma=1700 if metaread==520 & zsh>0.8374373
replace zpuma=1800 if metaread==520 & zsh>0.8659716
replace zpuma=1900 if metaread==520 & zsh>0.9091496
replace zpuma=2000 if metaread==520 & zsh>0.9447584
replace zpuma=101 if metaread==560 & zsh>0
replace zpuma=102 if metaread==560 & zsh>0.3275463
replace zpuma=200 if metaread==560 & zsh>0.7223002
replace zpuma=900 if metaread==600 & zsh>0
replace zpuma=2300 if metaread==600 & zsh>0.3235382
replace zpuma=2400 if metaread==600 & zsh>0.7510114
replace zpuma=5201 if metaread==640 & zsh>0
replace zpuma=5202 if metaread==640 & zsh>0.1019996
replace zpuma=5301 if metaread==640 & zsh>0.2037727
replace zpuma=5302 if metaread==640 & zsh>0.3550742
replace zpuma=5303 if metaread==640 & zsh>0.4931883
replace zpuma=5304 if metaread==640 & zsh>0.6254949
replace zpuma=5401 if metaread==640 & zsh>0.7733329
replace zpuma=5402 if metaread==640 & zsh>0.9099254
replace zpuma=3800 if metaread==680 & zsh>0
replace zpuma=3901 if metaread==680 & zsh>0.3870593
replace zpuma=3902 if metaread==680 & zsh>0.5922581
replace zpuma=3903 if metaread==680 & zsh>0.7935306
replace zpuma=400 if metaread==720 & zsh>0
replace zpuma=501 if metaread==720 & zsh>0.0626889
replace zpuma=502 if metaread==720 & zsh>0.110804
replace zpuma=503 if metaread==720 & zsh>0.1552506
replace zpuma=504 if metaread==720 & zsh>0.1994321
replace zpuma=505 if metaread==720 & zsh>0.2467361
replace zpuma=506 if metaread==720 & zsh>0.2887712
replace zpuma=507 if metaread==720 & zsh>0.3316676
replace zpuma=601 if metaread==720 & zsh>0.3741663
replace zpuma=602 if metaread==720 & zsh>0.4224577
replace zpuma=801 if metaread==720 & zsh>0.4627596
replace zpuma=802 if metaread==720 & zsh>0.4977147
replace zpuma=803 if metaread==720 & zsh>0.5408082
replace zpuma=804 if metaread==720 & zsh>0.583935
replace zpuma=805 if metaread==720 & zsh>0.6339572
replace zpuma=806 if metaread==720 & zsh>0.6668022
replace zpuma=901 if metaread==720 & zsh>0.7032357
replace zpuma=902 if metaread==720 & zsh>0.7478815
replace zpuma=1201 if metaread==720 & zsh>0.8055584
replace zpuma=1202 if metaread==720 & zsh>0.8549852
replace zpuma=1203 if metaread==720 & zsh>0.8975646
replace zpuma=1204 if metaread==720 & zsh>0.9496807
replace zpuma=1401 if metaread==760 & zsh>0
replace zpuma=1402 if metaread==760 & zsh>0.1587302
replace zpuma=1501 if metaread==760 & zsh>0.3451629
replace zpuma=1502 if metaread==760 & zsh>0.503467
replace zpuma=1600 if metaread==760 & zsh>0.7269835
replace zpuma=4200 if metaread==840 & zsh>0
replace zpuma=4300 if metaread==840 & zsh>0.3440558
replace zpuma=4400 if metaread==840 & zsh>0.6424962
replace zpuma=2200 if metaread==920 & zsh>0
replace zpuma=2300 if metaread==920 & zsh>0.5927306
replace zpuma=2601 if metaread==960 & zsh>0
replace zpuma=2602 if metaread==960 & zsh>0.5033826
replace zpuma=901 if metaread==1000 & zsh>0
replace zpuma=902 if metaread==1000 & zsh>0.1260003
replace zpuma=903 if metaread==1000 & zsh>0.263736
replace zpuma=904 if metaread==1000 & zsh>0.4694552
replace zpuma=905 if metaread==1000 & zsh>0.6876625
replace zpuma=1300 if metaread==1000 & zsh>0.8052043
replace zpuma=500 if metaread==1080 & zsh>0
replace zpuma=600 if metaread==1080 & zsh>0.2901423
replace zpuma=700 if metaread==1080 & zsh>0.7351592
replace zpuma=400 if metaread==1120 & zsh>0
replace zpuma=900 if metaread==1120 & zsh>0.0307639
replace zpuma=1000 if metaread==1120 & zsh>0.0621291
replace zpuma=1100 if metaread==1120 & zsh>0.0943259
replace zpuma=1200 if metaread==1120 & zsh>0.130402
replace zpuma=1300 if metaread==1120 & zsh>0.1625968
replace zpuma=1400 if metaread==1120 & zsh>0.1938116
replace zpuma=2400 if metaread==1120 & zsh>0.2249731
replace zpuma=2500 if metaread==1120 & zsh>0.2546616
replace zpuma=2600 if metaread==1120 & zsh>0.2865182
replace zpuma=2700 if metaread==1120 & zsh>0.3148516
replace zpuma=2800 if metaread==1120 & zsh>0.3736531
replace zpuma=2900 if metaread==1120 & zsh>0.405259
replace zpuma=3000 if metaread==1120 & zsh>0.4321704
replace zpuma=3100 if metaread==1120 & zsh>0.4650969
replace zpuma=3200 if metaread==1120 & zsh>0.5012363
replace zpuma=3301 if metaread==1120 & zsh>0.5366958
replace zpuma=3302 if metaread==1120 & zsh>0.5754458
replace zpuma=3303 if metaread==1120 & zsh>0.617855
replace zpuma=3304 if metaread==1120 & zsh>0.647416
replace zpuma=3305 if metaread==1120 & zsh>0.6758027
replace zpuma=3400 if metaread==1120 & zsh>0.7099861
replace zpuma=3500 if metaread==1120 & zsh>0.754905
replace zpuma=3600 if metaread==1120 & zsh>0.7848616
replace zpuma=3700 if metaread==1120 & zsh>0.815134
replace zpuma=3800 if metaread==1120 & zsh>0.8493803
replace zpuma=3900 if metaread==1120 & zsh>0.8832166
replace zpuma=4200 if metaread==1120 & zsh>0.9259195
replace zpuma=4600 if metaread==1120 & zsh>0.9574596
replace zpuma=700 if metaread==1121 & zsh>0
replace zpuma=800 if metaread==1121 & zsh>0.5275896
replace zpuma=500 if metaread==1122 & zsh>0
replace zpuma=600 if metaread==1122 & zsh>0.6531777
replace zpuma=1701 if metaread==1150 & zsh>0
replace zpuma=1702 if metaread==1150 & zsh>0.4891362
replace zpuma=2200 if metaread==1160 & zsh>0
replace zpuma=2400 if metaread==1160 & zsh>0.6100979
replace zpuma=4000 if metaread==1200 & zsh>0
replace zpuma=4100 if metaread==1200 & zsh>0.4240597
replace zpuma=6800 if metaread==1240 & zsh>0
replace zpuma=6900 if metaread==1240 & zsh>0.4201141
replace zpuma=1401 if metaread==1280 & zsh>0
replace zpuma=1402 if metaread==1280 & zsh>0.1397938
replace zpuma=1501 if metaread==1280 & zsh>0.2314938
replace zpuma=1502 if metaread==1280 & zsh>0.3269826
replace zpuma=1601 if metaread==1280 & zsh>0.4220088
replace zpuma=1602 if metaread==1280 & zsh>0.5197177
replace zpuma=1603 if metaread==1280 & zsh>0.6256506
replace zpuma=1604 if metaread==1280 & zsh>0.7139821
replace zpuma=1605 if metaread==1280 & zsh>0.8504308
replace zpuma=2201 if metaread==1320 & zsh>0
replace zpuma=2202 if metaread==1320 & zsh>0.3440816
replace zpuma=2203 if metaread==1320 & zsh>0.6906542
replace zpuma=2101 if metaread==1440 & zsh>0
replace zpuma=2102 if metaread==1440 & zsh>0.4132968
replace zpuma=2200 if metaread==1440 & zsh>0.7089472
replace zpuma=500 if metaread==1520 & zsh>0
replace zpuma=901 if metaread==1520 & zsh>0.1063024
replace zpuma=902 if metaread==1520 & zsh>0.1776845
replace zpuma=903 if metaread==1520 & zsh>0.2441583
replace zpuma=904 if metaread==1520 & zsh>0.3171352
replace zpuma=905 if metaread==1520 & zsh>0.3881949
replace zpuma=1000 if metaread==1520 & zsh>0.4781201
replace zpuma=1100 if metaread==1520 & zsh>0.5838667
replace zpuma=1200 if metaread==1520 & zsh>0.748041
replace zpuma=1300 if metaread==1520 & zsh>0.8309624
replace zpuma=1400 if metaread==1520 & zsh>0.9163259
replace zpuma=100 if metaread==1560 & zsh>0
replace zpuma=1800 if metaread==1560 & zsh>0.2952804
replace zpuma=1900 if metaread==1560 & zsh>0.641161
replace zpuma=3001 if metaread==1600 & zsh>0
replace zpuma=3002 if metaread==1600 & zsh>0.0143382
replace zpuma=3003 if metaread==1600 & zsh>0.0289949
replace zpuma=3004 if metaread==1600 & zsh>0.0476162
replace zpuma=3005 if metaread==1600 & zsh>0.0614054
replace zpuma=3006 if metaread==1600 & zsh>0.0760793
replace zpuma=3101 if metaread==1600 & zsh>0.0927446
replace zpuma=3102 if metaread==1600 & zsh>0.1083137
replace zpuma=3103 if metaread==1600 & zsh>0.1216166
replace zpuma=3104 if metaread==1600 & zsh>0.1397177
replace zpuma=3201 if metaread==1600 & zsh>0.159531
replace zpuma=3202 if metaread==1600 & zsh>0.1855548
replace zpuma=3203 if metaread==1600 & zsh>0.2008452
replace zpuma=3204 if metaread==1600 & zsh>0.2167251
replace zpuma=3205 if metaread==1600 & zsh>0.2355256
replace zpuma=3206 if metaread==1600 & zsh>0.2515811
replace zpuma=3301 if metaread==1600 & zsh>0.2785886
replace zpuma=3302 if metaread==1600 & zsh>0.2904007
replace zpuma=3303 if metaread==1600 & zsh>0.3063166
replace zpuma=3304 if metaread==1600 & zsh>0.3312955
replace zpuma=3305 if metaread==1600 & zsh>0.3449992
replace zpuma=3401 if metaread==1600 & zsh>0.3588099
replace zpuma=3402 if metaread==1600 & zsh>0.3761581
replace zpuma=3403 if metaread==1600 & zsh>0.4055116
replace zpuma=3404 if metaread==1600 & zsh>0.4394392
replace zpuma=3405 if metaread==1600 & zsh>0.4667622
replace zpuma=3406 if metaread==1600 & zsh>0.4961462
replace zpuma=3407 if metaread==1600 & zsh>0.5114458
replace zpuma=3408 if metaread==1600 & zsh>0.5332302
replace zpuma=3409 if metaread==1600 & zsh>0.5564548
replace zpuma=3410 if metaread==1600 & zsh>0.5746927
replace zpuma=3411 if metaread==1600 & zsh>0.5948172
replace zpuma=3412 if metaread==1600 & zsh>0.6153269
replace zpuma=3413 if metaread==1600 & zsh>0.6285595
replace zpuma=3414 if metaread==1600 & zsh>0.6486138
replace zpuma=3501 if metaread==1600 & zsh>0.6673556
replace zpuma=3502 if metaread==1600 & zsh>0.6931372
replace zpuma=3503 if metaread==1600 & zsh>0.7203799
replace zpuma=3504 if metaread==1600 & zsh>0.7396714
replace zpuma=3505 if metaread==1600 & zsh>0.7581769
replace zpuma=3506 if metaread==1600 & zsh>0.7734256
replace zpuma=3507 if metaread==1600 & zsh>0.7927058
replace zpuma=3508 if metaread==1600 & zsh>0.8037537
replace zpuma=3509 if metaread==1600 & zsh>0.8172914
replace zpuma=3510 if metaread==1600 & zsh>0.845866
replace zpuma=3511 if metaread==1600 & zsh>0.8663543
replace zpuma=3512 if metaread==1600 & zsh>0.8820332
replace zpuma=3513 if metaread==1600 & zsh>0.9005223
replace zpuma=3514 if metaread==1600 & zsh>0.9216376
replace zpuma=3515 if metaread==1600 & zsh>0.9351902
replace zpuma=3516 if metaread==1600 & zsh>0.9524196
replace zpuma=3517 if metaread==1600 & zsh>0.9675529
replace zpuma=3518 if metaread==1600 & zsh>0.9788175
replace zpuma=3519 if metaread==1600 & zsh>0.9889982
replace zpuma=100 if metaread==1602 & zsh>0
replace zpuma=201 if metaread==1602 & zsh>0.1412575
replace zpuma=202 if metaread==1602 & zsh>0.3170148
replace zpuma=203 if metaread==1602 & zsh>0.5775939
replace zpuma=300 if metaread==1602 & zsh>0.7531044
replace zpuma=2300 if metaread==1640 & zsh>0
replace zpuma=2400 if metaread==1640 & zsh>0.0691286
replace zpuma=4200 if metaread==1640 & zsh>0.1728458
replace zpuma=4401 if metaread==1640 & zsh>0.2821762
replace zpuma=4402 if metaread==1640 & zsh>0.3557703
replace zpuma=4403 if metaread==1640 & zsh>0.4550529
replace zpuma=4404 if metaread==1640 & zsh>0.5503701
replace zpuma=4501 if metaread==1640 & zsh>0.624026
replace zpuma=4502 if metaread==1640 & zsh>0.7022948
replace zpuma=4503 if metaread==1640 & zsh>0.7733597
replace zpuma=4601 if metaread==1640 & zsh>0.8478234
replace zpuma=4602 if metaread==1640 & zsh>0.9233254
replace zpuma=501 if metaread==1680 & zsh>0
replace zpuma=502 if metaread==1680 & zsh>0.0606827
replace zpuma=601 if metaread==1680 & zsh>0.1284282
replace zpuma=602 if metaread==1680 & zsh>0.1827619
replace zpuma=603 if metaread==1680 & zsh>0.2473085
replace zpuma=604 if metaread==1680 & zsh>0.2923413
replace zpuma=605 if metaread==1680 & zsh>0.345715
replace zpuma=606 if metaread==1680 & zsh>0.4029063
replace zpuma=607 if metaread==1680 & zsh>0.4463021
replace zpuma=608 if metaread==1680 & zsh>0.4969654
replace zpuma=609 if metaread==1680 & zsh>0.546508
replace zpuma=610 if metaread==1680 & zsh>0.594786
replace zpuma=611 if metaread==1680 & zsh>0.6438043
replace zpuma=612 if metaread==1680 & zsh>0.6900234
replace zpuma=701 if metaread==1680 & zsh>0.7358299
replace zpuma=702 if metaread==1680 & zsh>0.7853078
replace zpuma=703 if metaread==1680 & zsh>0.8349591
replace zpuma=800 if metaread==1680 & zsh>0.8851811
replace zpuma=1300 if metaread==1680 & zsh>0.9293427
replace zpuma=901 if metaread==1720 & zsh>0
replace zpuma=902 if metaread==1720 & zsh>0.1771076
replace zpuma=903 if metaread==1720 & zsh>0.390826
replace zpuma=904 if metaread==1720 & zsh>0.5975499
replace zpuma=905 if metaread==1720 & zsh>0.8158295
replace zpuma=1001 if metaread==1760 & zsh>0
replace zpuma=1002 if metaread==1760 & zsh>0.1909361
replace zpuma=1100 if metaread==1760 & zsh>0.406938
replace zpuma=1200 if metaread==1760 & zsh>0.6205807
replace zpuma=2700 if metaread==1840 & zsh>0
replace zpuma=3101 if metaread==1840 & zsh>0.0752605
replace zpuma=3102 if metaread==1840 & zsh>0.1460233
replace zpuma=3103 if metaread==1840 & zsh>0.2400316
replace zpuma=3104 if metaread==1840 & zsh>0.3072959
replace zpuma=3105 if metaread==1840 & zsh>0.4112294
replace zpuma=3106 if metaread==1840 & zsh>0.5072362
replace zpuma=3107 if metaread==1840 & zsh>0.5966022
replace zpuma=3108 if metaread==1840 & zsh>0.6605491
replace zpuma=3109 if metaread==1840 & zsh>0.748167
replace zpuma=3200 if metaread==1840 & zsh>0.8208871
replace zpuma=3300 if metaread==1840 & zsh>0.9186149
replace zpuma=6302 if metaread==1880 & zsh>0
replace zpuma=6400 if metaread==1880 & zsh>0.6252175
replace zpuma=2000 if metaread==1920 & zsh>0
replace zpuma=2101 if metaread==1920 & zsh>0.0314266
replace zpuma=2102 if metaread==1920 & zsh>0.0665902
replace zpuma=2103 if metaread==1920 & zsh>0.1161142
replace zpuma=2104 if metaread==1920 & zsh>0.1506925
replace zpuma=2201 if metaread==1920 & zsh>0.184245
replace zpuma=2202 if metaread==1920 & zsh>0.2604378
replace zpuma=2301 if metaread==1920 & zsh>0.3216292
replace zpuma=2302 if metaread==1920 & zsh>0.3494087
replace zpuma=2303 if metaread==1920 & zsh>0.392031
replace zpuma=2304 if metaread==1920 & zsh>0.4258132
replace zpuma=2305 if metaread==1920 & zsh>0.480578
replace zpuma=2306 if metaread==1920 & zsh>0.526432
replace zpuma=2307 if metaread==1920 & zsh>0.5843152
replace zpuma=2308 if metaread==1920 & zsh>0.6203316
replace zpuma=2309 if metaread==1920 & zsh>0.662796
replace zpuma=2310 if metaread==1920 & zsh>0.6967003
replace zpuma=2311 if metaread==1920 & zsh>0.7278461
replace zpuma=2312 if metaread==1920 & zsh>0.7889456
replace zpuma=2313 if metaread==1920 & zsh>0.8275149
replace zpuma=2314 if metaread==1920 & zsh>0.857637
replace zpuma=2315 if metaread==1920 & zsh>0.9031894
replace zpuma=2400 if metaread==1920 & zsh>0.9687022
replace zpuma=2501 if metaread==1921 & zsh>0
replace zpuma=2502 if metaread==1921 & zsh>0.0766113
replace zpuma=2503 if metaread==1921 & zsh>0.1436152
replace zpuma=2504 if metaread==1921 & zsh>0.2579668
replace zpuma=2505 if metaread==1921 & zsh>0.3532317
replace zpuma=2506 if metaread==1921 & zsh>0.4516584
replace zpuma=2507 if metaread==1921 & zsh>0.5654983
replace zpuma=2508 if metaread==1921 & zsh>0.6553441
replace zpuma=2509 if metaread==1921 & zsh>0.7180122
replace zpuma=2510 if metaread==1921 & zsh>0.8028743
replace zpuma=2511 if metaread==1921 & zsh>0.8691138
replace zpuma=2600 if metaread==1921 & zsh>0.9268044
replace zpuma=103 if metaread==1960 & zsh>0
replace zpuma=1900 if metaread==1960 & zsh>0.3978602
replace zpuma=3001 if metaread==2000 & zsh>0
replace zpuma=3002 if metaread==2000 & zsh>0.1131801
replace zpuma=3900 if metaread==2000 & zsh>0.2615667
replace zpuma=4000 if metaread==2000 & zsh>0.4164138
replace zpuma=4101 if metaread==2000 & zsh>0.5832828
replace zpuma=4102 if metaread==2000 & zsh>0.7099624
replace zpuma=4103 if metaread==2000 & zsh>0.8336393
replace zpuma=1901 if metaread==2020 & zsh>0
replace zpuma=1902 if metaread==2020 & zsh>0.4356141
replace zpuma=1903 if metaread==2020 & zsh>0.7189634
replace zpuma=805 if metaread==2080 & zsh>0
replace zpuma=806 if metaread==2080 & zsh>0.0550602
replace zpuma=807 if metaread==2080 & zsh>0.1110183
replace zpuma=808 if metaread==2080 & zsh>0.1594749
replace zpuma=809 if metaread==2080 & zsh>0.2163418
replace zpuma=810 if metaread==2080 & zsh>0.2685559
replace zpuma=811 if metaread==2080 & zsh>0.3246447
replace zpuma=812 if metaread==2080 & zsh>0.3759315
replace zpuma=813 if metaread==2080 & zsh>0.4313083
replace zpuma=814 if metaread==2080 & zsh>0.4800614
replace zpuma=815 if metaread==2080 & zsh>0.5415201
replace zpuma=816 if metaread==2080 & zsh>0.6194131
replace zpuma=817 if metaread==2080 & zsh>0.6718366
replace zpuma=818 if metaread==2080 & zsh>0.725823
replace zpuma=819 if metaread==2080 & zsh>0.7760662
replace zpuma=820 if metaread==2080 & zsh>0.8432984
replace zpuma=821 if metaread==2080 & zsh>0.8965402
replace zpuma=822 if metaread==2080 & zsh>0.9514102
replace zpuma=802 if metaread==2081 & zsh>0
replace zpuma=803 if metaread==2081 & zsh>0.4636092
replace zpuma=1400 if metaread==2120 & zsh>0
replace zpuma=1500 if metaread==2120 & zsh>0.5211749
replace zpuma=2501 if metaread==2160 & zsh>0
replace zpuma=2502 if metaread==2160 & zsh>0.0371287
replace zpuma=2503 if metaread==2160 & zsh>0.0838905
replace zpuma=2504 if metaread==2160 & zsh>0.1227529
replace zpuma=2505 if metaread==2160 & zsh>0.1621662
replace zpuma=2506 if metaread==2160 & zsh>0.1976347
replace zpuma=2507 if metaread==2160 & zsh>0.2391046
replace zpuma=2508 if metaread==2160 & zsh>0.2633855
replace zpuma=3100 if metaread==2160 & zsh>0.2898701
replace zpuma=3500 if metaread==2160 & zsh>0.3235183
replace zpuma=3600 if metaread==2160 & zsh>0.3813641
replace zpuma=3701 if metaread==2160 & zsh>0.4048134
replace zpuma=3702 if metaread==2160 & zsh>0.4287364
replace zpuma=3703 if metaread==2160 & zsh>0.4488982
replace zpuma=3704 if metaread==2160 & zsh>0.4696084
replace zpuma=3705 if metaread==2160 & zsh>0.4954284
replace zpuma=3706 if metaread==2160 & zsh>0.517131
replace zpuma=3707 if metaread==2160 & zsh>0.5388511
replace zpuma=3708 if metaread==2160 & zsh>0.5573291
replace zpuma=3801 if metaread==2160 & zsh>0.5854032
replace zpuma=3802 if metaread==2160 & zsh>0.6066996
replace zpuma=3803 if metaread==2160 & zsh>0.6495759
replace zpuma=3804 if metaread==2160 & zsh>0.689518
replace zpuma=3805 if metaread==2160 & zsh>0.7158621
replace zpuma=3806 if metaread==2160 & zsh>0.7451134
replace zpuma=3807 if metaread==2160 & zsh>0.7845218
replace zpuma=3900 if metaread==2160 & zsh>0.8136647
replace zpuma=4000 if metaread==2160 & zsh>0.8437228
replace zpuma=4101 if metaread==2160 & zsh>0.8750535
replace zpuma=4102 if metaread==2160 & zsh>0.9105874
replace zpuma=4103 if metaread==2160 & zsh>0.9478847
replace zpuma=4104 if metaread==2160 & zsh>0.9732863
replace zpuma=3201 if metaread==2281 & zsh>0
replace zpuma=3202 if metaread==2281 & zsh>0.4343806
replace zpuma=3501 if metaread==2310 & zsh>0
replace zpuma=3502 if metaread==2310 & zsh>0.1616189
replace zpuma=3503 if metaread==2310 & zsh>0.3292309
replace zpuma=3504 if metaread==2310 & zsh>0.4776301
replace zpuma=3505 if metaread==2310 & zsh>0.7475789
replace zpuma=100 if metaread==2360 & zsh>0
replace zpuma=200 if metaread==2360 & zsh>0.3572007
replace zpuma=701 if metaread==2400 & zsh>0
replace zpuma=702 if metaread==2400 & zsh>0.3151349
replace zpuma=3700 if metaread==2440 & zsh>0
replace zpuma=3800 if metaread==2440 & zsh>0.4861155
replace zpuma=3600 if metaread==2560 & zsh>0
replace zpuma=3700 if metaread==2560 & zsh>0.4214145
replace zpuma=100 if metaread==2580 & zsh>0
replace zpuma=200 if metaread==2580 & zsh>0.4734381
replace zpuma=2200 if metaread==2640 & zsh>0
replace zpuma=2301 if metaread==2640 & zsh>0.5062472
replace zpuma=102 if metaread==2670 & zsh>0
replace zpuma=200 if metaread==2670 & zsh>0.4763176
replace zpuma=3601 if metaread==2680 & zsh>0
replace zpuma=3602 if metaread==2680 & zsh>0.0848759
replace zpuma=3603 if metaread==2680 & zsh>0.1459121
replace zpuma=3604 if metaread==2680 & zsh>0.2109135
replace zpuma=3605 if metaread==2680 & zsh>0.2810207
replace zpuma=3606 if metaread==2680 & zsh>0.3753005
replace zpuma=3607 if metaread==2680 & zsh>0.4792465
replace zpuma=3608 if metaread==2680 & zsh>0.5416491
replace zpuma=3609 if metaread==2680 & zsh>0.6306528
replace zpuma=3610 if metaread==2680 & zsh>0.7019095
replace zpuma=3611 if metaread==2680 & zsh>0.7859788
replace zpuma=3612 if metaread==2680 & zsh>0.8646812
replace zpuma=3613 if metaread==2680 & zsh>0.9233786
replace zpuma=3800 if metaread==2700 & zsh>0
replace zpuma=3901 if metaread==2700 & zsh>0.2440581
replace zpuma=3902 if metaread==2700 & zsh>0.5100138
replace zpuma=3903 if metaread==2700 & zsh>0.7392428
replace zpuma=3300 if metaread==2710 & zsh>0
replace zpuma=3400 if metaread==2710 & zsh>0.6002974
replace zpuma=900 if metaread==2760 & zsh>0
replace zpuma=1000 if metaread==2760 & zsh>0.4421673
replace zpuma=1100 if metaread==2760 & zsh>0.7213189
replace zpuma=3200 if metaread==2840 & zsh>0
replace zpuma=3301 if metaread==2840 & zsh>0.1286645
replace zpuma=3302 if metaread==2840 & zsh>0.4158075
replace zpuma=3401 if metaread==2840 & zsh>0.5754626
replace zpuma=3402 if metaread==2840 & zsh>0.7360309
replace zpuma=1001 if metaread==2900 & zsh>0
replace zpuma=1002 if metaread==2900 & zsh>0.5153707
replace zpuma=4901 if metaread==2920 & zsh>0
replace zpuma=4902 if metaread==2920 & zsh>0.4480184
replace zpuma=1200 if metaread==3000 & zsh>0
replace zpuma=1300 if metaread==3000 & zsh>0.1616432
replace zpuma=1401 if metaread==3000 & zsh>0.3623488
replace zpuma=1402 if metaread==3000 & zsh>0.5013554
replace zpuma=1403 if metaread==3000 & zsh>0.6093032
replace zpuma=1500 if metaread==3000 & zsh>0.7482313
replace zpuma=200 if metaread==3080 & zsh>0
replace zpuma=300 if metaread==3080 & zsh>0.4443843
replace zpuma=1601 if metaread==3120 & zsh>0
replace zpuma=1602 if metaread==3120 & zsh>0.0867579
replace zpuma=1700 if metaread==3120 & zsh>0.1857112
replace zpuma=1800 if metaread==3120 & zsh>0.3428477
replace zpuma=1900 if metaread==3120 & zsh>0.4876985
replace zpuma=2000 if metaread==3120 & zsh>0.5853637
replace zpuma=2100 if metaread==3120 & zsh>0.6898383
replace zpuma=2200 if metaread==3120 & zsh>0.7772754
replace zpuma=2300 if metaread==3120 & zsh>0.8942227
replace zpuma=201 if metaread==3160 & zsh>0
replace zpuma=202 if metaread==3160 & zsh>0.2992706
replace zpuma=301 if metaread==3160 & zsh>0.4896468
replace zpuma=302 if metaread==3160 & zsh>0.6221935
replace zpuma=700 if metaread==3160 & zsh>0.8020127
replace zpuma=4301 if metaread==3200 & zsh>0
replace zpuma=4302 if metaread==3200 & zsh>0.4240918
replace zpuma=2900 if metaread==3240 & zsh>0
replace zpuma=3001 if metaread==3240 & zsh>0.2034543
replace zpuma=3002 if metaread==3240 & zsh>0.3932325
replace zpuma=3101 if metaread==3240 & zsh>0.5842131
replace zpuma=3102 if metaread==3240 & zsh>0.791896
replace zpuma=300 if metaread==3280 & zsh>0
replace zpuma=400 if metaread==3280 & zsh>0.1493221
replace zpuma=500 if metaread==3280 & zsh>0.3696691
replace zpuma=600 if metaread==3280 & zsh>0.5341617
replace zpuma=700 if metaread==3280 & zsh>0.6887341
replace zpuma=800 if metaread==3280 & zsh>0.8510188
replace zpuma=700 if metaread==3290 & zsh>0
replace zpuma=800 if metaread==3290 & zsh>0.4217198
replace zpuma=301 if metaread==3320 & zsh>0
replace zpuma=302 if metaread==3320 & zsh>0.1105358
replace zpuma=303 if metaread==3320 & zsh>0.2401963
replace zpuma=304 if metaread==3320 & zsh>0.3683551
replace zpuma=305 if metaread==3320 & zsh>0.5158947
replace zpuma=306 if metaread==3320 & zsh>0.6465067
replace zpuma=307 if metaread==3320 & zsh>0.7960532
replace zpuma=4501 if metaread==3360 & zsh>0
replace zpuma=4502 if metaread==3360 & zsh>0.0308577
replace zpuma=4503 if metaread==3360 & zsh>0.0714101
replace zpuma=4601 if metaread==3360 & zsh>0.0990374
replace zpuma=4602 if metaread==3360 & zsh>0.1261047
replace zpuma=4603 if metaread==3360 & zsh>0.1559956
replace zpuma=4604 if metaread==3360 & zsh>0.1862779
replace zpuma=4605 if metaread==3360 & zsh>0.2240053
replace zpuma=4606 if metaread==3360 & zsh>0.2525018
replace zpuma=4607 if metaread==3360 & zsh>0.2821183
replace zpuma=4608 if metaread==3360 & zsh>0.3182354
replace zpuma=4609 if metaread==3360 & zsh>0.3416506
replace zpuma=4610 if metaread==3360 & zsh>0.3728299
replace zpuma=4611 if metaread==3360 & zsh>0.4039864
replace zpuma=4612 if metaread==3360 & zsh>0.432825
replace zpuma=4613 if metaread==3360 & zsh>0.4643331
replace zpuma=4614 if metaread==3360 & zsh>0.4968764
replace zpuma=4615 if metaread==3360 & zsh>0.5309252
replace zpuma=4616 if metaread==3360 & zsh>0.563141
replace zpuma=4617 if metaread==3360 & zsh>0.5956891
replace zpuma=4618 if metaread==3360 & zsh>0.6326535
replace zpuma=4619 if metaread==3360 & zsh>0.674441
replace zpuma=4620 if metaread==3360 & zsh>0.7126099
replace zpuma=4621 if metaread==3360 & zsh>0.7393312
replace zpuma=4622 if metaread==3360 & zsh>0.7792966
replace zpuma=4623 if metaread==3360 & zsh>0.8134373
replace zpuma=4624 if metaread==3360 & zsh>0.8505062
replace zpuma=4625 if metaread==3360 & zsh>0.8838609
replace zpuma=4701 if metaread==3360 & zsh>0.9157672
replace zpuma=4702 if metaread==3360 & zsh>0.9565054
replace zpuma=4801 if metaread==3361 & zsh>0
replace zpuma=4802 if metaread==3361 & zsh>0.5253463
replace zpuma=200 if metaread==3440 & zsh>0
replace zpuma=300 if metaread==3440 & zsh>0.4708771
replace zpuma=1901 if metaread==3480 & zsh>0
replace zpuma=1902 if metaread==3480 & zsh>0.0718437
replace zpuma=2001 if metaread==3480 & zsh>0.1418531
replace zpuma=2002 if metaread==3480 & zsh>0.2185349
replace zpuma=2301 if metaread==3480 & zsh>0.2825518
replace zpuma=2302 if metaread==3480 & zsh>0.3510497
replace zpuma=2303 if metaread==3480 & zsh>0.4203899
replace zpuma=2304 if metaread==3480 & zsh>0.4875794
replace zpuma=2305 if metaread==3480 & zsh>0.5723627
replace zpuma=2306 if metaread==3480 & zsh>0.6724864
replace zpuma=2307 if metaread==3480 & zsh>0.7491264
replace zpuma=2400 if metaread==3480 & zsh>0.8223127
replace zpuma=2500 if metaread==3480 & zsh>0.8938026
replace zpuma=1300 if metaread==3560 & zsh>0
replace zpuma=1400 if metaread==3560 & zsh>0.2738548
replace zpuma=1500 if metaread==3560 & zsh>0.6702494
replace zpuma=1101 if metaread==3590 & zsh>0
replace zpuma=1102 if metaread==3590 & zsh>0.0949435
replace zpuma=1103 if metaread==3590 & zsh>0.2046371
replace zpuma=1104 if metaread==3590 & zsh>0.3045698
replace zpuma=1105 if metaread==3590 & zsh>0.4169509
replace zpuma=1106 if metaread==3590 & zsh>0.5242494
replace zpuma=1107 if metaread==3590 & zsh>0.6212421
replace zpuma=1200 if metaread==3590 & zsh>0.7602971
replace zpuma=1300 if metaread==3590 & zsh>0.8767314
replace zpuma=801 if metaread==3660 & zsh>0
replace zpuma=802 if metaread==3660 & zsh>0.333149
replace zpuma=1000 if metaread==3660 & zsh>0.6424626
replace zpuma=2501 if metaread==3680 & zsh>0
replace zpuma=2502 if metaread==3680 & zsh>0.5554879
replace zpuma=2601 if metaread==3720 & zsh>0
replace zpuma=2602 if metaread==3720 & zsh>0.3028281
replace zpuma=2800 if metaread==3720 & zsh>0.7165608
replace zpuma=500 if metaread==3760 & zsh>0
replace zpuma=601 if metaread==3760 & zsh>0.0874365
replace zpuma=602 if metaread==3760 & zsh>0.1726501
replace zpuma=603 if metaread==3760 & zsh>0.2372294
replace zpuma=604 if metaread==3760 & zsh>0.3059477
replace zpuma=800 if metaread==3760 & zsh>0.3690105
replace zpuma=901 if metaread==3760 & zsh>0.465609
replace zpuma=902 if metaread==3760 & zsh>0.5719948
replace zpuma=1001 if metaread==3760 & zsh>0.6766925
replace zpuma=1002 if metaread==3760 & zsh>0.7509431
replace zpuma=1003 if metaread==3760 & zsh>0.8087056
replace zpuma=1004 if metaread==3760 & zsh>0.87689
replace zpuma=1100 if metaread==3760 & zsh>0.9342963
replace zpuma=3701 if metaread==3810 & zsh>0
replace zpuma=3702 if metaread==3810 & zsh>0.4082856
replace zpuma=3703 if metaread==3810 & zsh>0.7303595
replace zpuma=1301 if metaread==3840 & zsh>0
replace zpuma=1302 if metaread==3840 & zsh>0.337109
replace zpuma=1400 if metaread==3840 & zsh>0.5031362
replace zpuma=1500 if metaread==3840 & zsh>0.8171151
replace zpuma=1000 if metaread==3880 & zsh>0
replace zpuma=1100 if metaread==3880 & zsh>0.4671837
replace zpuma=2401 if metaread==3980 & zsh>0
replace zpuma=2402 if metaread==3980 & zsh>0.3186223
replace zpuma=2403 if metaread==3980 & zsh>0.5764576
replace zpuma=2404 if metaread==3980 & zsh>0.7698063
replace zpuma=3301 if metaread==4000 & zsh>0
replace zpuma=3302 if metaread==4000 & zsh>0.3492152
replace zpuma=3303 if metaread==4000 & zsh>0.6559213
replace zpuma=1800 if metaread==4040 & zsh>0
replace zpuma=1900 if metaread==4040 & zsh>0.2546005
replace zpuma=2000 if metaread==4040 & zsh>0.6576876
replace zpuma=501 if metaread==4120 & zsh>0
replace zpuma=502 if metaread==4120 & zsh>0.0829213
replace zpuma=503 if metaread==4120 & zsh>0.1479373
replace zpuma=504 if metaread==4120 & zsh>0.2281671
replace zpuma=505 if metaread==4120 & zsh>0.3387881
replace zpuma=506 if metaread==4120 & zsh>0.4130165
replace zpuma=507 if metaread==4120 & zsh>0.4928094
replace zpuma=508 if metaread==4120 & zsh>0.5732457
replace zpuma=509 if metaread==4120 & zsh>0.6727352
replace zpuma=510 if metaread==4120 & zsh>0.7673888
replace zpuma=511 if metaread==4120 & zsh>0.8649317
replace zpuma=1901 if metaread==4280 & zsh>0
replace zpuma=1902 if metaread==4280 & zsh>0.4271137
replace zpuma=801 if metaread==4360 & zsh>0
replace zpuma=802 if metaread==4360 & zsh>0.5669383
replace zpuma=900 if metaread==4400 & zsh>0
replace zpuma=1000 if metaread==4400 & zsh>0.3173282
replace zpuma=1100 if metaread==4400 & zsh>0.6128311
replace zpuma=4500 if metaread==4480 & zsh>0
replace zpuma=4600 if metaread==4480 & zsh>0.0108163
replace zpuma=4700 if metaread==4480 & zsh>0.0214486
replace zpuma=4800 if metaread==4480 & zsh>0.0384361
replace zpuma=4900 if metaread==4480 & zsh>0.05005
replace zpuma=5000 if metaread==4480 & zsh>0.0705839
replace zpuma=5100 if metaread==4480 & zsh>0.0856963
replace zpuma=5200 if metaread==4480 & zsh>0.0966296
replace zpuma=5300 if metaread==4480 & zsh>0.1112913
replace zpuma=5401 if metaread==4480 & zsh>0.1225388
replace zpuma=5402 if metaread==4480 & zsh>0.1420043
replace zpuma=5403 if metaread==4480 & zsh>0.1558169
replace zpuma=5404 if metaread==4480 & zsh>0.1672224
replace zpuma=5405 if metaread==4480 & zsh>0.1811541
replace zpuma=5406 if metaread==4480 & zsh>0.193869
replace zpuma=5407 if metaread==4480 & zsh>0.2087921
replace zpuma=5408 if metaread==4480 & zsh>0.226647
replace zpuma=5409 if metaread==4480 & zsh>0.2444359
replace zpuma=5410 if metaread==4480 & zsh>0.2622704
replace zpuma=5411 if metaread==4480 & zsh>0.2828309
replace zpuma=5412 if metaread==4480 & zsh>0.3090436
replace zpuma=5413 if metaread==4480 & zsh>0.3288147
replace zpuma=5414 if metaread==4480 & zsh>0.3568431
replace zpuma=5415 if metaread==4480 & zsh>0.3703603
replace zpuma=5416 if metaread==4480 & zsh>0.3914575
replace zpuma=5417 if metaread==4480 & zsh>0.4149024
replace zpuma=5418 if metaread==4480 & zsh>0.429471
replace zpuma=5419 if metaread==4480 & zsh>0.4400167
replace zpuma=5420 if metaread==4480 & zsh>0.4547268
replace zpuma=5421 if metaread==4480 & zsh>0.4698539
replace zpuma=5422 if metaread==4480 & zsh>0.4833287
replace zpuma=5423 if metaread==4480 & zsh>0.4936669
replace zpuma=5424 if metaread==4480 & zsh>0.5036478
replace zpuma=5500 if metaread==4480 & zsh>0.5226966
replace zpuma=5600 if metaread==4480 & zsh>0.533488
replace zpuma=5701 if metaread==4480 & zsh>0.5490491
replace zpuma=5702 if metaread==4480 & zsh>0.5625917
replace zpuma=5703 if metaread==4480 & zsh>0.5808766
replace zpuma=5800 if metaread==4480 & zsh>0.596674
replace zpuma=5900 if metaread==4480 & zsh>0.6076107
replace zpuma=6000 if metaread==4480 & zsh>0.6187023
replace zpuma=6101 if metaread==4480 & zsh>0.6288994
replace zpuma=6102 if metaread==4480 & zsh>0.6416581
replace zpuma=6103 if metaread==4480 & zsh>0.6570005
replace zpuma=6104 if metaread==4480 & zsh>0.6689945
replace zpuma=6105 if metaread==4480 & zsh>0.6869942
replace zpuma=6106 if metaread==4480 & zsh>0.7028103
replace zpuma=6107 if metaread==4480 & zsh>0.7210112
replace zpuma=6108 if metaread==4480 & zsh>0.7326173
replace zpuma=6109 if metaread==4480 & zsh>0.7445416
replace zpuma=6110 if metaread==4480 & zsh>0.7570968
replace zpuma=6111 if metaread==4480 & zsh>0.7739482
replace zpuma=6112 if metaread==4480 & zsh>0.785763
replace zpuma=6113 if metaread==4480 & zsh>0.7980719
replace zpuma=6114 if metaread==4480 & zsh>0.8117237
replace zpuma=6115 if metaread==4480 & zsh>0.8247421
replace zpuma=6116 if metaread==4480 & zsh>0.8402075
replace zpuma=6117 if metaread==4480 & zsh>0.8525077
replace zpuma=6118 if metaread==4480 & zsh>0.8656398
replace zpuma=6119 if metaread==4480 & zsh>0.8830807
replace zpuma=6120 if metaread==4480 & zsh>0.8954174
replace zpuma=6121 if metaread==4480 & zsh>0.9065467
replace zpuma=6122 if metaread==4480 & zsh>0.9173808
replace zpuma=6123 if metaread==4480 & zsh>0.9348446
replace zpuma=6124 if metaread==4480 & zsh>0.9485322
replace zpuma=6125 if metaread==4480 & zsh>0.9597681
replace zpuma=6126 if metaread==4480 & zsh>0.9879295
replace zpuma=6801 if metaread==4482 & zsh>0
replace zpuma=6802 if metaread==4482 & zsh>0.0622148
replace zpuma=6900 if metaread==4482 & zsh>0.1114252
replace zpuma=7000 if metaread==4482 & zsh>0.1545015
replace zpuma=7100 if metaread==4482 & zsh>0.199373
replace zpuma=7200 if metaread==4482 & zsh>0.2536578
replace zpuma=7300 if metaread==4482 & zsh>0.3252601
replace zpuma=7400 if metaread==4482 & zsh>0.3793246
replace zpuma=7501 if metaread==4482 & zsh>0.4248082
replace zpuma=7502 if metaread==4482 & zsh>0.4739486
replace zpuma=7601 if metaread==4482 & zsh>0.5311056
replace zpuma=7602 if metaread==4482 & zsh>0.5941056
replace zpuma=7603 if metaread==4482 & zsh>0.6537454
replace zpuma=7604 if metaread==4482 & zsh>0.7232455
replace zpuma=7605 if metaread==4482 & zsh>0.7951523
replace zpuma=7606 if metaread==4482 & zsh>0.8574979
replace zpuma=7607 if metaread==4482 & zsh>0.9272699
replace zpuma=1701 if metaread==4520 & zsh>0
replace zpuma=1702 if metaread==4520 & zsh>0.1414427
replace zpuma=1703 if metaread==4520 & zsh>0.2594242
replace zpuma=1704 if metaread==4520 & zsh>0.483632
replace zpuma=1705 if metaread==4520 & zsh>0.6011391
replace zpuma=3300 if metaread==4520 & zsh>0.7547033
replace zpuma=3400 if metaread==4520 & zsh>0.8837455
replace zpuma=501 if metaread==4600 & zsh>0
replace zpuma=502 if metaread==4600 & zsh>0.4641403
replace zpuma=2700 if metaread==4680 & zsh>0
replace zpuma=2800 if metaread==4680 & zsh>0.4962109
replace zpuma=1100 if metaread==4720 & zsh>0
replace zpuma=1200 if metaread==4720 & zsh>0.5150418
replace zpuma=6600 if metaread==4880 & zsh>0
replace zpuma=6701 if metaread==4880 & zsh>0.2114892
replace zpuma=6702 if metaread==4880 & zsh>0.4010066
replace zpuma=6703 if metaread==4880 & zsh>0.5803972
replace zpuma=6704 if metaread==4880 & zsh>0.7925302
replace zpuma=2101 if metaread==4900 & zsh>0
replace zpuma=2102 if metaread==4900 & zsh>0.32462
replace zpuma=2103 if metaread==4900 & zsh>0.5783982
replace zpuma=2104 if metaread==4900 & zsh>0.793079
replace zpuma=100 if metaread==4920 & zsh>0
replace zpuma=3101 if metaread==4920 & zsh>0.1109676
replace zpuma=3102 if metaread==4920 & zsh>0.2597492
replace zpuma=3103 if metaread==4920 & zsh>0.3949398
replace zpuma=3104 if metaread==4920 & zsh>0.5166777
replace zpuma=3105 if metaread==4920 & zsh>0.6441435
replace zpuma=3201 if metaread==4920 & zsh>0.7384711
replace zpuma=3202 if metaread==4920 & zsh>0.8888371
replace zpuma=4001 if metaread==5000 & zsh>0
replace zpuma=4002 if metaread==5000 & zsh>0.0435414
replace zpuma=4003 if metaread==5000 & zsh>0.0873712
replace zpuma=4004 if metaread==5000 & zsh>0.1385667
replace zpuma=4005 if metaread==5000 & zsh>0.1903407
replace zpuma=4006 if metaread==5000 & zsh>0.2444591
replace zpuma=4007 if metaread==5000 & zsh>0.2847212
replace zpuma=4008 if metaread==5000 & zsh>0.3400372
replace zpuma=4009 if metaread==5000 & zsh>0.4016175
replace zpuma=4010 if metaread==5000 & zsh>0.4439183
replace zpuma=4011 if metaread==5000 & zsh>0.4908028
replace zpuma=4012 if metaread==5000 & zsh>0.5648108
replace zpuma=4013 if metaread==5000 & zsh>0.6325642
replace zpuma=4014 if metaread==5000 & zsh>0.6851288
replace zpuma=4015 if metaread==5000 & zsh>0.7391238
replace zpuma=4016 if metaread==5000 & zsh>0.8025223
replace zpuma=4017 if metaread==5000 & zsh>0.8521051
replace zpuma=4018 if metaread==5000 & zsh>0.9071236
replace zpuma=4019 if metaread==5000 & zsh>0.9514953
replace zpuma=2001 if metaread==5080 & zsh>0
replace zpuma=2002 if metaread==5080 & zsh>0.1142475
replace zpuma=2003 if metaread==5080 & zsh>0.1987362
replace zpuma=2004 if metaread==5080 & zsh>0.2747459
replace zpuma=2101 if metaread==5080 & zsh>0.3750158
replace zpuma=2102 if metaread==5080 & zsh>0.4516312
replace zpuma=2201 if metaread==5080 & zsh>0.6098399
replace zpuma=2202 if metaread==5080 & zsh>0.7112248
replace zpuma=2203 if metaread==5080 & zsh>0.7818874
replace zpuma=2300 if metaread==5080 & zsh>0.8633922
replace zpuma=900 if metaread==5120 & zsh>0
replace zpuma=1001 if metaread==5120 & zsh>0.0783103
replace zpuma=1002 if metaread==5120 & zsh>0.1343496
replace zpuma=1100 if metaread==5120 & zsh>0.1814418
replace zpuma=1201 if metaread==5120 & zsh>0.2371
replace zpuma=1202 if metaread==5120 & zsh>0.2873701
replace zpuma=1203 if metaread==5120 & zsh>0.3258671
replace zpuma=1301 if metaread==5120 & zsh>0.3617145
replace zpuma=1302 if metaread==5120 & zsh>0.4070072
replace zpuma=1303 if metaread==5120 & zsh>0.4574239
replace zpuma=1401 if metaread==5120 & zsh>0.5010398
replace zpuma=1402 if metaread==5120 & zsh>0.5384003
replace zpuma=1403 if metaread==5120 & zsh>0.6060991
replace zpuma=1404 if metaread==5120 & zsh>0.6434751
replace zpuma=1405 if metaread==5120 & zsh>0.6820411
replace zpuma=1406 if metaread==5120 & zsh>0.719975
replace zpuma=1501 if metaread==5120 & zsh>0.7633484
replace zpuma=1502 if metaread==5120 & zsh>0.8192885
replace zpuma=1601 if metaread==5120 & zsh>0.8559592
replace zpuma=1602 if metaread==5120 & zsh>0.8914934
replace zpuma=1700 if metaread==5120 & zsh>0.9328318
replace zpuma=2200 if metaread==5160 & zsh>0
replace zpuma=2300 if metaread==5160 & zsh>0.3680435
replace zpuma=2400 if metaread==5160 & zsh>0.7352827
replace zpuma=2500 if metaread==5170 & zsh>0
replace zpuma=2601 if metaread==5170 & zsh>0.4243543
replace zpuma=2602 if metaread==5170 & zsh>0.7280025
replace zpuma=1101 if metaread==5190 & zsh>0
replace zpuma=1102 if metaread==5190 & zsh>0.105819
replace zpuma=1103 if metaread==5190 & zsh>0.241394
replace zpuma=1104 if metaread==5190 & zsh>0.3528827
replace zpuma=1105 if metaread==5190 & zsh>0.4701232
replace zpuma=1201 if metaread==5190 & zsh>0.5742072
replace zpuma=1202 if metaread==5190 & zsh>0.7287912
replace zpuma=1203 if metaread==5190 & zsh>0.8598141
replace zpuma=1900 if metaread==5240 & zsh>0
replace zpuma=2000 if metaread==5240 & zsh>0.6005007
replace zpuma=3701 if metaread==5340 & zsh>0
replace zpuma=3702 if metaread==5340 & zsh>0.4765286
replace zpuma=400 if metaread==5360 & zsh>0
replace zpuma=501 if metaread==5360 & zsh>0.1047959
replace zpuma=502 if metaread==5360 & zsh>0.1956097
replace zpuma=2201 if metaread==5360 & zsh>0.2788779
replace zpuma=2202 if metaread==5360 & zsh>0.3823518
replace zpuma=2203 if metaread==5360 & zsh>0.4603267
replace zpuma=2204 if metaread==5360 & zsh>0.5698119
replace zpuma=2205 if metaread==5360 & zsh>0.6629241
replace zpuma=2300 if metaread==5360 & zsh>0.746558
replace zpuma=2400 if metaread==5360 & zsh>0.8986253
replace zpuma=1500 if metaread==5480 & zsh>0
replace zpuma=1800 if metaread==5480 & zsh>0.3470773
replace zpuma=2000 if metaread==5480 & zsh>0.6630037
replace zpuma=1801 if metaread==5560 & zsh>0
replace zpuma=1802 if metaread==5560 & zsh>0.0813899
replace zpuma=1803 if metaread==5560 & zsh>0.1706978
replace zpuma=1804 if metaread==5560 & zsh>0.2719372
replace zpuma=1901 if metaread==5560 & zsh>0.3722396
replace zpuma=1902 if metaread==5560 & zsh>0.4993464
replace zpuma=1903 if metaread==5560 & zsh>0.5933068
replace zpuma=1904 if metaread==5560 & zsh>0.674961
replace zpuma=1905 if metaread==5560 & zsh>0.7555798
replace zpuma=2002 if metaread==5560 & zsh>0.8324049
replace zpuma=2400 if metaread==5560 & zsh>0.9150097
replace zpuma=3400 if metaread==5600 & zsh>0
replace zpuma=3501 if metaread==5600 & zsh>0.0211847
replace zpuma=3502 if metaread==5600 & zsh>0.0356925
replace zpuma=3503 if metaread==5600 & zsh>0.0518215
replace zpuma=3504 if metaread==5600 & zsh>0.0706431
replace zpuma=3505 if metaread==5600 & zsh>0.0836957
replace zpuma=3506 if metaread==5600 & zsh>0.1044157
replace zpuma=3601 if metaread==5600 & zsh>0.1163406
replace zpuma=3602 if metaread==5600 & zsh>0.1319769
replace zpuma=3701 if metaread==5600 & zsh>0.1484842
replace zpuma=3702 if metaread==5600 & zsh>0.1594806
replace zpuma=3703 if metaread==5600 & zsh>0.1738137
replace zpuma=3704 if metaread==5600 & zsh>0.1850636
replace zpuma=3705 if metaread==5600 & zsh>0.1972365
replace zpuma=3706 if metaread==5600 & zsh>0.2088319
replace zpuma=3707 if metaread==5600 & zsh>0.2212476
replace zpuma=3708 if metaread==5600 & zsh>0.2330767
replace zpuma=3709 if metaread==5600 & zsh>0.2444628
replace zpuma=3710 if metaread==5600 & zsh>0.2609272
replace zpuma=3801 if metaread==5600 & zsh>0.2733293
replace zpuma=3802 if metaread==5600 & zsh>0.2974525
replace zpuma=3803 if metaread==5600 & zsh>0.3116944
replace zpuma=3804 if metaread==5600 & zsh>0.3212154
replace zpuma=3805 if metaread==5600 & zsh>0.3314173
replace zpuma=3806 if metaread==5600 & zsh>0.3629951
replace zpuma=3807 if metaread==5600 & zsh>0.3903218
replace zpuma=3808 if metaread==5600 & zsh>0.408857
replace zpuma=3809 if metaread==5600 & zsh>0.4306489
replace zpuma=3810 if metaread==5600 & zsh>0.4500807
replace zpuma=3901 if metaread==5600 & zsh>0.4695354
replace zpuma=3902 if metaread==5600 & zsh>0.4863468
replace zpuma=3903 if metaread==5600 & zsh>0.5000626
replace zpuma=4001 if metaread==5600 & zsh>0.5167744
replace zpuma=4002 if metaread==5600 & zsh>0.5302739
replace zpuma=4003 if metaread==5600 & zsh>0.5406611
replace zpuma=4004 if metaread==5600 & zsh>0.5507067
replace zpuma=4005 if metaread==5600 & zsh>0.5642821
replace zpuma=4006 if metaread==5600 & zsh>0.5781459
replace zpuma=4007 if metaread==5600 & zsh>0.590107
replace zpuma=4008 if metaread==5600 & zsh>0.599255
replace zpuma=4009 if metaread==5600 & zsh>0.6116974
replace zpuma=4010 if metaread==5600 & zsh>0.6317352
replace zpuma=4011 if metaread==5600 & zsh>0.6470909
replace zpuma=4012 if metaread==5600 & zsh>0.6588069
replace zpuma=4013 if metaread==5600 & zsh>0.6732416
replace zpuma=4014 if metaread==5600 & zsh>0.6863332
replace zpuma=4015 if metaread==5600 & zsh>0.6997175
replace zpuma=4016 if metaread==5600 & zsh>0.7156617
replace zpuma=4017 if metaread==5600 & zsh>0.730361
replace zpuma=4018 if metaread==5600 & zsh>0.7476901
replace zpuma=4101 if metaread==5600 & zsh>0.7579612
replace zpuma=4102 if metaread==5600 & zsh>0.7799303
replace zpuma=4103 if metaread==5600 & zsh>0.7996666
replace zpuma=4104 if metaread==5600 & zsh>0.827103
replace zpuma=4105 if metaread==5600 & zsh>0.8400275
replace zpuma=4106 if metaread==5600 & zsh>0.8605359
replace zpuma=4107 if metaread==5600 & zsh>0.8771579
replace zpuma=4108 if metaread==5600 & zsh>0.8935015
replace zpuma=4109 if metaread==5600 & zsh>0.9065788
replace zpuma=4110 if metaread==5600 & zsh>0.9228862
replace zpuma=4111 if metaread==5600 & zsh>0.9406981
replace zpuma=4112 if metaread==5600 & zsh>0.9553408
replace zpuma=4113 if metaread==5600 & zsh>0.9769748
replace zpuma=4114 if metaread==5600 & zsh>0.9904188
replace zpuma=4201 if metaread==5601 & zsh>0
replace zpuma=4202 if metaread==5601 & zsh>0.0404702
replace zpuma=4203 if metaread==5601 & zsh>0.0817756
replace zpuma=4204 if metaread==5601 & zsh>0.1195796
replace zpuma=4205 if metaread==5601 & zsh>0.1598126
replace zpuma=4206 if metaread==5601 & zsh>0.1983179
replace zpuma=4207 if metaread==5601 & zsh>0.2491947
replace zpuma=4208 if metaread==5601 & zsh>0.2911409
replace zpuma=4209 if metaread==5601 & zsh>0.3272797
replace zpuma=4210 if metaread==5601 & zsh>0.3656718
replace zpuma=4211 if metaread==5601 & zsh>0.4035483
replace zpuma=4212 if metaread==5601 & zsh>0.4421782
replace zpuma=4301 if metaread==5601 & zsh>0.4819537
replace zpuma=4302 if metaread==5601 & zsh>0.5529363
replace zpuma=4303 if metaread==5601 & zsh>0.5942071
replace zpuma=4304 if metaread==5601 & zsh>0.6364031
replace zpuma=4305 if metaread==5601 & zsh>0.6817868
replace zpuma=4306 if metaread==5601 & zsh>0.724122
replace zpuma=4307 if metaread==5601 & zsh>0.7665184
replace zpuma=4308 if metaread==5601 & zsh>0.8057657
replace zpuma=4309 if metaread==5601 & zsh>0.8446946
replace zpuma=4310 if metaread==5601 & zsh>0.8862452
replace zpuma=4311 if metaread==5601 & zsh>0.9246233
replace zpuma=4312 if metaread==5601 & zsh>0.9607694
replace zpuma=301 if metaread==5602 & zsh>0
replace zpuma=302 if metaread==5602 & zsh>0.0877111
replace zpuma=303 if metaread==5602 & zsh>0.1742209
replace zpuma=304 if metaread==5602 & zsh>0.2707099
replace zpuma=305 if metaread==5602 & zsh>0.4156215
replace zpuma=306 if metaread==5602 & zsh>0.5428368
replace zpuma=400 if metaread==5602 & zsh>0.6512876
replace zpuma=501 if metaread==5602 & zsh>0.7502506
replace zpuma=502 if metaread==5602 & zsh>0.8524932
replace zpuma=601 if metaread==5603 & zsh>0
replace zpuma=602 if metaread==5603 & zsh>0.2111851
replace zpuma=701 if metaread==5603 & zsh>0.3796147
replace zpuma=702 if metaread==5603 & zsh>0.5933606
replace zpuma=703 if metaread==5603 & zsh>0.8104868
replace zpuma=800 if metaread==5604 & zsh>0
replace zpuma=901 if metaread==5604 & zsh>0.1047999
replace zpuma=902 if metaread==5604 & zsh>0.2267906
replace zpuma=903 if metaread==5604 & zsh>0.3697398
replace zpuma=904 if metaread==5604 & zsh>0.4938242
replace zpuma=905 if metaread==5604 & zsh>0.6060472
replace zpuma=1001 if metaread==5604 & zsh>0.7420939
replace zpuma=1002 if metaread==5604 & zsh>0.865701
replace zpuma=1301 if metaread==5605 & zsh>0
replace zpuma=1302 if metaread==5605 & zsh>0.0666342
replace zpuma=1401 if metaread==5605 & zsh>0.1159138
replace zpuma=1402 if metaread==5605 & zsh>0.1752058
replace zpuma=1403 if metaread==5605 & zsh>0.2314949
replace zpuma=1404 if metaread==5605 & zsh>0.3005397
replace zpuma=1501 if metaread==5605 & zsh>0.3731963
replace zpuma=1502 if metaread==5605 & zsh>0.4356816
replace zpuma=1503 if metaread==5605 & zsh>0.4979405
replace zpuma=1504 if metaread==5605 & zsh>0.5588287
replace zpuma=1600 if metaread==5605 & zsh>0.6190211
replace zpuma=1700 if metaread==5605 & zsh>0.6921357
replace zpuma=1800 if metaread==5605 & zsh>0.7426385
replace zpuma=1901 if metaread==5605 & zsh>0.7987916
replace zpuma=1902 if metaread==5605 & zsh>0.8795451
replace zpuma=1903 if metaread==5605 & zsh>0.9395707
replace zpuma=3301 if metaread==5660 & zsh>0
replace zpuma=3302 if metaread==5660 & zsh>0.3406663
replace zpuma=3303 if metaread==5660 & zsh>0.6759883
replace zpuma=2000 if metaread==5720 & zsh>0
replace zpuma=2100 if metaread==5720 & zsh>0.116019
replace zpuma=2200 if metaread==5720 & zsh>0.2072043
replace zpuma=2801 if metaread==5720 & zsh>0.3236312
replace zpuma=2802 if metaread==5720 & zsh>0.4014769
replace zpuma=2803 if metaread==5720 & zsh>0.5040508
replace zpuma=2900 if metaread==5720 & zsh>0.6050301
replace zpuma=3000 if metaread==5720 & zsh>0.7450294
replace zpuma=3100 if metaread==5720 & zsh>0.8759902
replace zpuma=1501 if metaread==5790 & zsh>0
replace zpuma=1502 if metaread==5790 & zsh>0.4369562
replace zpuma=3200 if metaread==5800 & zsh>0
replace zpuma=3300 if metaread==5800 & zsh>0.4966612
replace zpuma=1301 if metaread==5880 & zsh>0
replace zpuma=1302 if metaread==5880 & zsh>0.3194235
replace zpuma=1400 if metaread==5880 & zsh>0.5682374
replace zpuma=1500 if metaread==5880 & zsh>0.8110757
replace zpuma=1201 if metaread==5910 & zsh>0
replace zpuma=1202 if metaread==5910 & zsh>0.4745895
replace zpuma=702 if metaread==5920 & zsh>0
replace zpuma=901 if metaread==5920 & zsh>0.1914515
replace zpuma=902 if metaread==5920 & zsh>0.4239946
replace zpuma=903 if metaread==5920 & zsh>0.6358338
replace zpuma=904 if metaread==5920 & zsh>0.8125207
replace zpuma=1801 if metaread==5960 & zsh>0
replace zpuma=1802 if metaread==5960 & zsh>0.0578297
replace zpuma=2001 if metaread==5960 & zsh>0.1086122
replace zpuma=2002 if metaread==5960 & zsh>0.2073881
replace zpuma=2003 if metaread==5960 & zsh>0.2707039
replace zpuma=2201 if metaread==5960 & zsh>0.3353157
replace zpuma=2202 if metaread==5960 & zsh>0.3953153
replace zpuma=2203 if metaread==5960 & zsh>0.5094974
replace zpuma=2204 if metaread==5960 & zsh>0.5900571
replace zpuma=2205 if metaread==5960 & zsh>0.6754301
replace zpuma=2206 if metaread==5960 & zsh>0.7379433
replace zpuma=2207 if metaread==5960 & zsh>0.8336388
replace zpuma=2300 if metaread==5960 & zsh>0.8984361
replace zpuma=101 if metaread==6080 & zsh>0
replace zpuma=102 if metaread==6080 & zsh>0.4436141
replace zpuma=200 if metaread==6080 & zsh>0.7062061
replace zpuma=1700 if metaread==6120 & zsh>0
replace zpuma=1800 if metaread==6120 & zsh>0.3177138
replace zpuma=1900 if metaread==6120 & zsh>0.6256149
replace zpuma=2001 if metaread==6160 & zsh>0
replace zpuma=2002 if metaread==6160 & zsh>0.0233291
replace zpuma=2003 if metaread==6160 & zsh>0.0541741
replace zpuma=2101 if metaread==6160 & zsh>0.0835946
replace zpuma=2102 if metaread==6160 & zsh>0.1061712
replace zpuma=2103 if metaread==6160 & zsh>0.1311212
replace zpuma=2104 if metaread==6160 & zsh>0.1600901
replace zpuma=2201 if metaread==6160 & zsh>0.1836469
replace zpuma=2202 if metaread==6160 & zsh>0.2081392
replace zpuma=3901 if metaread==6160 & zsh>0.2484694
replace zpuma=3902 if metaread==6160 & zsh>0.2888334
replace zpuma=3903 if metaread==6160 & zsh>0.3196856
replace zpuma=3904 if metaread==6160 & zsh>0.3470631
replace zpuma=4001 if metaread==6160 & zsh>0.3729819
replace zpuma=4002 if metaread==6160 & zsh>0.3963412
replace zpuma=4003 if metaread==6160 & zsh>0.4315254
replace zpuma=4004 if metaread==6160 & zsh>0.4645247
replace zpuma=4005 if metaread==6160 & zsh>0.4877479
replace zpuma=4006 if metaread==6160 & zsh>0.5108082
replace zpuma=4101 if metaread==6160 & zsh>0.5339926
replace zpuma=4102 if metaread==6160 & zsh>0.5637828
replace zpuma=4103 if metaread==6160 & zsh>0.5859979
replace zpuma=4104 if metaread==6160 & zsh>0.6057209
replace zpuma=4105 if metaread==6160 & zsh>0.6345189
replace zpuma=4106 if metaread==6160 & zsh>0.658617
replace zpuma=4107 if metaread==6160 & zsh>0.6864954
replace zpuma=4108 if metaread==6160 & zsh>0.7024032
replace zpuma=4109 if metaread==6160 & zsh>0.7240124
replace zpuma=4110 if metaread==6160 & zsh>0.7506352
replace zpuma=4111 if metaread==6160 & zsh>0.7698004
replace zpuma=4201 if metaread==6160 & zsh>0.7971274
replace zpuma=4202 if metaread==6160 & zsh>0.8249235
replace zpuma=4203 if metaread==6160 & zsh>0.8570122
replace zpuma=4204 if metaread==6160 & zsh>0.8859167
replace zpuma=4301 if metaread==6160 & zsh>0.9077573
replace zpuma=4302 if metaread==6160 & zsh>0.9371268
replace zpuma=4303 if metaread==6160 & zsh>0.97397
replace zpuma=101 if metaread==6200 & zsh>0
replace zpuma=102 if metaread==6200 & zsh>0.0189449
replace zpuma=103 if metaread==6200 & zsh>0.0530594
replace zpuma=104 if metaread==6200 & zsh>0.1138481
replace zpuma=105 if metaread==6200 & zsh>0.1814902
replace zpuma=106 if metaread==6200 & zsh>0.2362744
replace zpuma=107 if metaread==6200 & zsh>0.270751
replace zpuma=108 if metaread==6200 & zsh>0.3246611
replace zpuma=109 if metaread==6200 & zsh>0.3641919
replace zpuma=110 if metaread==6200 & zsh>0.3992766
replace zpuma=111 if metaread==6200 & zsh>0.429866
replace zpuma=112 if metaread==6200 & zsh>0.4803612
replace zpuma=113 if metaread==6200 & zsh>0.5370811
replace zpuma=114 if metaread==6200 & zsh>0.5717404
replace zpuma=115 if metaread==6200 & zsh>0.6276835
replace zpuma=116 if metaread==6200 & zsh>0.6816199
replace zpuma=117 if metaread==6200 & zsh>0.746567
replace zpuma=118 if metaread==6200 & zsh>0.7894706
replace zpuma=119 if metaread==6200 & zsh>0.8298532
replace zpuma=120 if metaread==6200 & zsh>0.8654922
replace zpuma=121 if metaread==6200 & zsh>0.9294234
replace zpuma=122 if metaread==6200 & zsh>0.9715244
replace zpuma=1701 if metaread==6280 & zsh>0
replace zpuma=1702 if metaread==6280 & zsh>0.0494231
replace zpuma=1703 if metaread==6280 & zsh>0.0976401
replace zpuma=1801 if metaread==6280 & zsh>0.1509054
replace zpuma=1802 if metaread==6280 & zsh>0.2160315
replace zpuma=1803 if metaread==6280 & zsh>0.2644075
replace zpuma=1804 if metaread==6280 & zsh>0.3222385
replace zpuma=1805 if metaread==6280 & zsh>0.4079075
replace zpuma=1806 if metaread==6280 & zsh>0.4628023
replace zpuma=1807 if metaread==6280 & zsh>0.5089388
replace zpuma=1900 if metaread==6280 & zsh>0.5702932
replace zpuma=2002 if metaread==6280 & zsh>0.6480872
replace zpuma=2101 if metaread==6280 & zsh>0.7191717
replace zpuma=2102 if metaread==6280 & zsh>0.7653171
replace zpuma=2103 if metaread==6280 & zsh>0.8323769
replace zpuma=2201 if metaread==6280 & zsh>0.8797901
replace zpuma=2300 if metaread==6280 & zsh>0.9420351
replace zpuma=200 if metaread==6400 & zsh>0
replace zpuma=300 if metaread==6400 & zsh>0.4375682
replace zpuma=1301 if metaread==6440 & zsh>0
replace zpuma=1302 if metaread==6440 & zsh>0.0563719
replace zpuma=1303 if metaread==6440 & zsh>0.1107247
replace zpuma=1304 if metaread==6440 & zsh>0.1685075
replace zpuma=1305 if metaread==6440 & zsh>0.2336688
replace zpuma=1306 if metaread==6440 & zsh>0.3005684
replace zpuma=1307 if metaread==6440 & zsh>0.3646995
replace zpuma=1308 if metaread==6440 & zsh>0.42067
replace zpuma=1309 if metaread==6440 & zsh>0.4884847
replace zpuma=1310 if metaread==6440 & zsh>0.548965
replace zpuma=1311 if metaread==6440 & zsh>0.6119462
replace zpuma=1312 if metaread==6440 & zsh>0.6932368
replace zpuma=1313 if metaread==6440 & zsh>0.7571077
replace zpuma=2101 if metaread==6440 & zsh>0.8171155
replace zpuma=2102 if metaread==6440 & zsh>0.8719464
replace zpuma=2200 if metaread==6440 & zsh>0.9243456
replace zpuma=100 if metaread==6480 & zsh>0
replace zpuma=200 if metaread==6480 & zsh>0.166505
replace zpuma=300 if metaread==6480 & zsh>0.3174396
replace zpuma=400 if metaread==6480 & zsh>0.4310882
replace zpuma=600 if metaread==6480 & zsh>0.6021293
replace zpuma=4300 if metaread==6480 & zsh>0.7701065
replace zpuma=4400 if metaread==6480 & zsh>0.8844722
replace zpuma=601 if metaread==6520 & zsh>0
replace zpuma=602 if metaread==6520 & zsh>0.3589076
replace zpuma=603 if metaread==6520 & zsh>0.7077942
replace zpuma=2601 if metaread==6640 & zsh>0
replace zpuma=2602 if metaread==6640 & zsh>0.1271792
replace zpuma=2701 if metaread==6640 & zsh>0.2448803
replace zpuma=2702 if metaread==6640 & zsh>0.3366375
replace zpuma=2703 if metaread==6640 & zsh>0.4347558
replace zpuma=2801 if metaread==6640 & zsh>0.5286523
replace zpuma=2802 if metaread==6640 & zsh>0.6329257
replace zpuma=2900 if metaread==6640 & zsh>0.7214912
replace zpuma=3000 if metaread==6640 & zsh>0.8697149
replace zpuma=3401 if metaread==6680 & zsh>0
replace zpuma=3402 if metaread==6680 & zsh>0.4939923
replace zpuma=100 if metaread==6720 & zsh>0
replace zpuma=200 if metaread==6720 & zsh>0.4637965
replace zpuma=1200 if metaread==6760 & zsh>0
replace zpuma=1300 if metaread==6760 & zsh>0.1958463
replace zpuma=1400 if metaread==6760 & zsh>0.4680677
replace zpuma=1500 if metaread==6760 & zsh>0.7337417
replace zpuma=2300 if metaread==6760 & zsh>0.8659626
replace zpuma=4000 if metaread==6780 & zsh>0
replace zpuma=4100 if metaread==6780 & zsh>0.0505908
replace zpuma=4200 if metaread==6780 & zsh>0.0962173
replace zpuma=4300 if metaread==6780 & zsh>0.1339028
replace zpuma=4401 if metaread==6780 & zsh>0.1867561
replace zpuma=4402 if metaread==6780 & zsh>0.2383807
replace zpuma=4403 if metaread==6780 & zsh>0.273764
replace zpuma=4404 if metaread==6780 & zsh>0.3333396
replace zpuma=4405 if metaread==6780 & zsh>0.3948916
replace zpuma=4406 if metaread==6780 & zsh>0.4378295
replace zpuma=4407 if metaread==6780 & zsh>0.4909031
replace zpuma=7700 if metaread==6780 & zsh>0.5209565
replace zpuma=7801 if metaread==6780 & zsh>0.563364
replace zpuma=7802 if metaread==6780 & zsh>0.5997729
replace zpuma=7900 if metaread==6780 & zsh>0.6476028
replace zpuma=8001 if metaread==6780 & zsh>0.6904112
replace zpuma=8002 if metaread==6780 & zsh>0.7791386
replace zpuma=8003 if metaread==6780 & zsh>0.8367719
replace zpuma=8004 if metaread==6780 & zsh>0.888309
replace zpuma=8005 if metaread==6780 & zsh>0.9514698
replace zpuma=1600 if metaread==6800 & zsh>0
replace zpuma=1700 if metaread==6800 & zsh>0.5037834
replace zpuma=901 if metaread==6840 & zsh>0
replace zpuma=902 if metaread==6840 & zsh>0.0930695
replace zpuma=1001 if metaread==6840 & zsh>0.1964683
replace zpuma=1002 if metaread==6840 & zsh>0.3104711
replace zpuma=1003 if metaread==6840 & zsh>0.4319308
replace zpuma=1004 if metaread==6840 & zsh>0.5307347
replace zpuma=1005 if metaread==6840 & zsh>0.7012006
replace zpuma=1100 if metaread==6840 & zsh>0.8028017
replace zpuma=1300 if metaread==6840 & zsh>0.902431
replace zpuma=2800 if metaread==6880 & zsh>0
replace zpuma=2900 if metaread==6880 & zsh>0.4461692
replace zpuma=1401 if metaread==6920 & zsh>0
replace zpuma=1402 if metaread==6920 & zsh>0.0777908
replace zpuma=1403 if metaread==6920 & zsh>0.1417008
replace zpuma=1501 if metaread==6920 & zsh>0.2324841
replace zpuma=1502 if metaread==6920 & zsh>0.3098277
replace zpuma=1503 if metaread==6920 & zsh>0.3802745
replace zpuma=1504 if metaread==6920 & zsh>0.4734385
replace zpuma=1505 if metaread==6920 & zsh>0.582522
replace zpuma=1506 if metaread==6920 & zsh>0.6730588
replace zpuma=1601 if metaread==6920 & zsh>0.7334884
replace zpuma=1602 if metaread==6920 & zsh>0.8227901
replace zpuma=1700 if metaread==6920 & zsh>0.8964466
replace zpuma=800 if metaread==6960 & zsh>0
replace zpuma=900 if metaread==6960 & zsh>0.4870349
replace zpuma=1101 if metaread==7040 & zsh>0
replace zpuma=1102 if metaread==7040 & zsh>0.0625061
replace zpuma=1201 if metaread==7040 & zsh>0.1140638
replace zpuma=1202 if metaread==7040 & zsh>0.1591826
replace zpuma=1500 if metaread==7040 & zsh>0.2218104
replace zpuma=1601 if metaread==7040 & zsh>0.2834056
replace zpuma=1602 if metaread==7040 & zsh>0.3404501
replace zpuma=1701 if metaread==7040 & zsh>0.3994553
replace zpuma=1702 if metaread==7040 & zsh>0.4459646
replace zpuma=1703 if metaread==7040 & zsh>0.4969861
replace zpuma=1704 if metaread==7040 & zsh>0.5420667
replace zpuma=1705 if metaread==7040 & zsh>0.587016
replace zpuma=1706 if metaread==7040 & zsh>0.6367395
replace zpuma=1707 if metaread==7040 & zsh>0.6835995
replace zpuma=1708 if metaread==7040 & zsh>0.7466411
replace zpuma=1801 if metaread==7040 & zsh>0.7980121
replace zpuma=1802 if metaread==7040 & zsh>0.8310795
replace zpuma=1803 if metaread==7040 & zsh>0.8702595
replace zpuma=1900 if metaread==7040 & zsh>0.9203974
replace zpuma=1101 if metaread==7080 & zsh>0
replace zpuma=1102 if metaread==7080 & zsh>0.6071422
replace zpuma=2900 if metaread==7120 & zsh>0
replace zpuma=3001 if metaread==7120 & zsh>0.5201219
replace zpuma=200 if metaread==7160 & zsh>0
replace zpuma=301 if metaread==7160 & zsh>0.1347604
replace zpuma=302 if metaread==7160 & zsh>0.2229834
replace zpuma=501 if metaread==7160 & zsh>0.2993484
replace zpuma=502 if metaread==7160 & zsh>0.4478352
replace zpuma=503 if metaread==7160 & zsh>0.5217986
replace zpuma=504 if metaread==7160 & zsh>0.6228386
replace zpuma=505 if metaread==7160 & zsh>0.7369617
replace zpuma=506 if metaread==7160 & zsh>0.8249429
replace zpuma=507 if metaread==7160 & zsh>0.9022522
replace zpuma=5500 if metaread==7240 & zsh>0
replace zpuma=5601 if metaread==7240 & zsh>0.1129149
replace zpuma=5602 if metaread==7240 & zsh>0.1712692
replace zpuma=5603 if metaread==7240 & zsh>0.2347207
replace zpuma=5604 if metaread==7240 & zsh>0.3005212
replace zpuma=5605 if metaread==7240 & zsh>0.3761753
replace zpuma=5606 if metaread==7240 & zsh>0.4628426
replace zpuma=5607 if metaread==7240 & zsh>0.5732313
replace zpuma=5608 if metaread==7240 & zsh>0.6977727
replace zpuma=5609 if metaread==7240 & zsh>0.7746174
replace zpuma=5610 if metaread==7240 & zsh>0.8521969
replace zpuma=5611 if metaread==7240 & zsh>0.9168724
replace zpuma=8101 if metaread==7320 & zsh>0
replace zpuma=8102 if metaread==7320 & zsh>0.0913681
replace zpuma=8103 if metaread==7320 & zsh>0.1326876
replace zpuma=8104 if metaread==7320 & zsh>0.186279
replace zpuma=8105 if metaread==7320 & zsh>0.2447921
replace zpuma=8106 if metaread==7320 & zsh>0.30696
replace zpuma=8107 if metaread==7320 & zsh>0.34373
replace zpuma=8108 if metaread==7320 & zsh>0.4036562
replace zpuma=8109 if metaread==7320 & zsh>0.486622
replace zpuma=8110 if metaread==7320 & zsh>0.5726513
replace zpuma=8111 if metaread==7320 & zsh>0.6303284
replace zpuma=8112 if metaread==7320 & zsh>0.6829462
replace zpuma=8113 if metaread==7320 & zsh>0.7430413
replace zpuma=8114 if metaread==7320 & zsh>0.8117034
replace zpuma=8115 if metaread==7320 & zsh>0.8732313
replace zpuma=8116 if metaread==7320 & zsh>0.9460946
replace zpuma=1201 if metaread==7360 & zsh>0
replace zpuma=1202 if metaread==7360 & zsh>0.0562771
replace zpuma=2201 if metaread==7360 & zsh>0.1438876
replace zpuma=2202 if metaread==7360 & zsh>0.2288638
replace zpuma=2203 if metaread==7360 & zsh>0.2991543
replace zpuma=2204 if metaread==7360 & zsh>0.363554
replace zpuma=2205 if metaread==7360 & zsh>0.4373458
replace zpuma=2206 if metaread==7360 & zsh>0.4986495
replace zpuma=2207 if metaread==7360 & zsh>0.5577417
replace zpuma=2301 if metaread==7360 & zsh>0.6075152
replace zpuma=2302 if metaread==7360 & zsh>0.6855537
replace zpuma=2303 if metaread==7360 & zsh>0.7478477
replace zpuma=2304 if metaread==7360 & zsh>0.8016509
replace zpuma=2305 if metaread==7360 & zsh>0.869552
replace zpuma=2306 if metaread==7360 & zsh>0.9409023
replace zpuma=2101 if metaread==7361 & zsh>0
replace zpuma=2102 if metaread==7361 & zsh>0.0466794
replace zpuma=2103 if metaread==7361 & zsh>0.0959152
replace zpuma=2104 if metaread==7361 & zsh>0.1428859
replace zpuma=2105 if metaread==7361 & zsh>0.2061187
replace zpuma=2106 if metaread==7361 & zsh>0.2542319
replace zpuma=2107 if metaread==7361 & zsh>0.3049689
replace zpuma=2108 if metaread==7361 & zsh>0.3462409
replace zpuma=2401 if metaread==7361 & zsh>0.390284
replace zpuma=2402 if metaread==7361 & zsh>0.4480758
replace zpuma=2403 if metaread==7361 & zsh>0.5125719
replace zpuma=2404 if metaread==7361 & zsh>0.5687902
replace zpuma=2405 if metaread==7361 & zsh>0.6138608
replace zpuma=2406 if metaread==7361 & zsh>0.6791554
replace zpuma=2407 if metaread==7361 & zsh>0.7335388
replace zpuma=2408 if metaread==7361 & zsh>0.7890815
replace zpuma=2409 if metaread==7361 & zsh>0.8428881
replace zpuma=2410 if metaread==7361 & zsh>0.9225749
replace zpuma=1000 if metaread==7362 & zsh>0
replace zpuma=1301 if metaread==7362 & zsh>0.2521025
replace zpuma=1302 if metaread==7362 & zsh>0.5375008
replace zpuma=1303 if metaread==7362 & zsh>0.760896
replace zpuma=2701 if metaread==7400 & zsh>0
replace zpuma=2702 if metaread==7400 & zsh>0.1139104
replace zpuma=2703 if metaread==7400 & zsh>0.203326
replace zpuma=2704 if metaread==7400 & zsh>0.2739968
replace zpuma=2705 if metaread==7400 & zsh>0.3472797
replace zpuma=2706 if metaread==7400 & zsh>0.4066834
replace zpuma=2707 if metaread==7400 & zsh>0.4625935
replace zpuma=2708 if metaread==7400 & zsh>0.5325272
replace zpuma=2709 if metaread==7400 & zsh>0.6088439
replace zpuma=2710 if metaread==7400 & zsh>0.6802508
replace zpuma=2711 if metaread==7400 & zsh>0.7653836
replace zpuma=2712 if metaread==7400 & zsh>0.8285266
replace zpuma=2713 if metaread==7400 & zsh>0.8901749
replace zpuma=2714 if metaread==7400 & zsh>0.9473979
replace zpuma=3701 if metaread==7460 & zsh>0
replace zpuma=3702 if metaread==7460 & zsh>0.448495
replace zpuma=6701 if metaread==7470 & zsh>0
replace zpuma=6702 if metaread==7470 & zsh>0.4323773
replace zpuma=2801 if metaread==7480 & zsh>0
replace zpuma=2802 if metaread==7480 & zsh>0.5082303
replace zpuma=1101 if metaread==7500 & zsh>0
replace zpuma=1102 if metaread==7500 & zsh>0.3783579
replace zpuma=1103 if metaread==7500 & zsh>0.6364081
replace zpuma=2801 if metaread==7510 & zsh>0
replace zpuma=2802 if metaread==7510 & zsh>0.2332571
replace zpuma=2901 if metaread==7510 & zsh>0.4481606
replace zpuma=2902 if metaread==7510 & zsh>0.7686378
replace zpuma=3400 if metaread==7520 & zsh>0
replace zpuma=3500 if metaread==7520 & zsh>0.533587
replace zpuma=801 if metaread==7560 & zsh>0
replace zpuma=802 if metaread==7560 & zsh>0.1765498
replace zpuma=901 if metaread==7560 & zsh>0.3839377
replace zpuma=902 if metaread==7560 & zsh>0.5491868
replace zpuma=903 if metaread==7560 & zsh>0.7322384
replace zpuma=1001 if metaread==7600 & zsh>0
replace zpuma=1002 if metaread==7600 & zsh>0.0537284
replace zpuma=1003 if metaread==7600 & zsh>0.1176733
replace zpuma=1004 if metaread==7600 & zsh>0.1607437
replace zpuma=1005 if metaread==7600 & zsh>0.2022257
replace zpuma=1801 if metaread==7600 & zsh>0.2446378
replace zpuma=1802 if metaread==7600 & zsh>0.3045284
replace zpuma=1803 if metaread==7600 & zsh>0.3568773
replace zpuma=1804 if metaread==7600 & zsh>0.4097261
replace zpuma=1805 if metaread==7600 & zsh>0.4542208
replace zpuma=1900 if metaread==7600 & zsh>0.5062606
replace zpuma=2001 if metaread==7600 & zsh>0.5550059
replace zpuma=2002 if metaread==7600 & zsh>0.605355
replace zpuma=2003 if metaread==7600 & zsh>0.667373
replace zpuma=2004 if metaread==7600 & zsh>0.7099177
replace zpuma=2005 if metaread==7600 & zsh>0.7586265
replace zpuma=2006 if metaread==7600 & zsh>0.8067459
replace zpuma=2007 if metaread==7600 & zsh>0.8547033
replace zpuma=2008 if metaread==7600 & zsh>0.9028507
replace zpuma=2009 if metaread==7600 & zsh>0.9495611
replace zpuma=101 if metaread==7680 & zsh>0
replace zpuma=102 if metaread==7680 & zsh>0.3428478
replace zpuma=200 if metaread==7680 & zsh>0.6415759
replace zpuma=500 if metaread==7800 & zsh>0
replace zpuma=600 if metaread==7800 & zsh>0.3786281
replace zpuma=500 if metaread==7840 & zsh>0
replace zpuma=601 if metaread==7840 & zsh>0.4731815
replace zpuma=602 if metaread==7840 & zsh>0.7177126
replace zpuma=2400 if metaread==7920 & zsh>0
replace zpuma=2500 if metaread==7920 & zsh>0.5251942
replace zpuma=1600 if metaread==8000 & zsh>0
replace zpuma=1700 if metaread==8000 & zsh>0.2053467
replace zpuma=1800 if metaread==8000 & zsh>0.376551
replace zpuma=1900 if metaread==8000 & zsh>0.575507
replace zpuma=2000 if metaread==8000 & zsh>0.8047649
replace zpuma=2300 if metaread==8040 & zsh>0
replace zpuma=2500 if metaread==8040 & zsh>0.6492581
replace zpuma=1901 if metaread==8120 & zsh>0
replace zpuma=1902 if metaread==8120 & zsh>0.1836137
replace zpuma=2001 if metaread==8120 & zsh>0.403434
replace zpuma=2002 if metaread==8120 & zsh>0.6863252
replace zpuma=600 if metaread==8160 & zsh>0
replace zpuma=700 if metaread==8160 & zsh>0.165075
replace zpuma=801 if metaread==8160 & zsh>0.3643341
replace zpuma=802 if metaread==8160 & zsh>0.5063394
replace zpuma=803 if metaread==8160 & zsh>0.683183
replace zpuma=804 if metaread==8160 & zsh>0.8617733
replace zpuma=1300 if metaread==8200 & zsh>0
replace zpuma=1401 if metaread==8200 & zsh>0.2819054
replace zpuma=1402 if metaread==8200 & zsh>0.5048734
replace zpuma=1403 if metaread==8200 & zsh>0.7230802
replace zpuma=1404 if metaread==8200 & zsh>0.8645232
replace zpuma=701 if metaread==8240 & zsh>0
replace zpuma=702 if metaread==8240 & zsh>0.6570795
replace zpuma=1700 if metaread==8280 & zsh>0
replace zpuma=2501 if metaread==8280 & zsh>0.0441261
replace zpuma=2502 if metaread==8280 & zsh>0.0890349
replace zpuma=2503 if metaread==8280 & zsh>0.1310384
replace zpuma=2601 if metaread==8280 & zsh>0.1699285
replace zpuma=2602 if metaread==8280 & zsh>0.2196892
replace zpuma=2603 if metaread==8280 & zsh>0.2616159
replace zpuma=2604 if metaread==8280 & zsh>0.3106093
replace zpuma=2605 if metaread==8280 & zsh>0.3554095
replace zpuma=2606 if metaread==8280 & zsh>0.4033906
replace zpuma=2607 if metaread==8280 & zsh>0.4509821
replace zpuma=2608 if metaread==8280 & zsh>0.5040609
replace zpuma=2701 if metaread==8280 & zsh>0.5595199
replace zpuma=2702 if metaread==8280 & zsh>0.6160037
replace zpuma=2703 if metaread==8280 & zsh>0.6833029
replace zpuma=2704 if metaread==8280 & zsh>0.7377813
replace zpuma=2705 if metaread==8280 & zsh>0.8069999
replace zpuma=2706 if metaread==8280 & zsh>0.8575022
replace zpuma=2707 if metaread==8280 & zsh>0.9007819
replace zpuma=2708 if metaread==8280 & zsh>0.959455
replace zpuma=201 if metaread==8400 & zsh>0
replace zpuma=202 if metaread==8400 & zsh>0.2581376
replace zpuma=301 if metaread==8400 & zsh>0.5097628
replace zpuma=302 if metaread==8400 & zsh>0.6821797
replace zpuma=303 if metaread==8400 & zsh>0.8573424
replace zpuma=2301 if metaread==8480 & zsh>0
replace zpuma=2302 if metaread==8480 & zsh>0.5441142
replace zpuma=201 if metaread==8520 & zsh>0
replace zpuma=202 if metaread==8520 & zsh>0.156444
replace zpuma=203 if metaread==8520 & zsh>0.2840102
replace zpuma=204 if metaread==8520 & zsh>0.4117796
replace zpuma=205 if metaread==8520 & zsh>0.5856392
replace zpuma=206 if metaread==8520 & zsh>0.7090222
replace zpuma=207 if metaread==8520 & zsh>0.8336383
replace zpuma=1000 if metaread==8560 & zsh>0
replace zpuma=1100 if metaread==8560 & zsh>0.1822663
replace zpuma=1200 if metaread==8560 & zsh>0.7475778
replace zpuma=401 if metaread==8680 & zsh>0
replace zpuma=402 if metaread==8680 & zsh>0.5260797
replace zpuma=6200 if metaread==8730 & zsh>0
replace zpuma=6300 if metaread==8730 & zsh>0.2171439
replace zpuma=6400 if metaread==8730 & zsh>0.3518244
replace zpuma=6500 if metaread==8730 & zsh>0.5131663
replace zpuma=6601 if metaread==8730 & zsh>0.6785633
replace zpuma=6602 if metaread==8730 & zsh>0.8354604
replace zpuma=3501 if metaread==8780 & zsh>0
replace zpuma=3502 if metaread==8780 & zsh>0.2959788
replace zpuma=3503 if metaread==8780 & zsh>0.6008001
replace zpuma=3801 if metaread==8800 & zsh>0
replace zpuma=3802 if metaread==8800 & zsh>0.4736972
replace zpuma=100 if metaread==8840 & zsh>0
replace zpuma=101 if metaread==8840 & zsh>0.0426584
replace zpuma=102 if metaread==8840 & zsh>0.0672762
replace zpuma=103 if metaread==8840 & zsh>0.0888464
replace zpuma=104 if metaread==8840 & zsh>0.110467
replace zpuma=105 if metaread==8840 & zsh>0.1319896
replace zpuma=200 if metaread==8840 & zsh>0.1566864
replace zpuma=300 if metaread==8840 & zsh>0.1859822
replace zpuma=301 if metaread==8840 & zsh>0.2299922
replace zpuma=302 if metaread==8840 & zsh>0.275626
replace zpuma=303 if metaread==8840 & zsh>0.3098927
replace zpuma=304 if metaread==8840 & zsh>0.3522286
replace zpuma=305 if metaread==8840 & zsh>0.3856785
replace zpuma=501 if metaread==8840 & zsh>0.4462732
replace zpuma=502 if metaread==8840 & zsh>0.4787761
replace zpuma=600 if metaread==8840 & zsh>0.5144595
replace zpuma=800 if metaread==8840 & zsh>0.5753617
replace zpuma=900 if metaread==8840 & zsh>0.5962713
replace zpuma=1001 if metaread==8840 & zsh>0.6270953
replace zpuma=1002 if metaread==8840 & zsh>0.6493348
replace zpuma=1003 if metaread==8840 & zsh>0.6756795
replace zpuma=1004 if metaread==8840 & zsh>0.7092873
replace zpuma=1005 if metaread==8840 & zsh>0.7443672
replace zpuma=1006 if metaread==8840 & zsh>0.769903
replace zpuma=1007 if metaread==8840 & zsh>0.7933993
replace zpuma=1101 if metaread==8840 & zsh>0.8171368
replace zpuma=1102 if metaread==8840 & zsh>0.8417942
replace zpuma=1103 if metaread==8840 & zsh>0.8677937
replace zpuma=1104 if metaread==8840 & zsh>0.8888141
replace zpuma=1105 if metaread==8840 & zsh>0.9064092
replace zpuma=1106 if metaread==8840 & zsh>0.9342186
replace zpuma=1107 if metaread==8840 & zsh>0.9580659
replace zpuma=1600 if metaread==8840 & zsh>0.9771733
replace zpuma=3501 if metaread==8960 & zsh>0
replace zpuma=3502 if metaread==8960 & zsh>0.1684314
replace zpuma=3503 if metaread==8960 & zsh>0.2836692
replace zpuma=3504 if metaread==8960 & zsh>0.3834049
replace zpuma=3505 if metaread==8960 & zsh>0.4809144
replace zpuma=3506 if metaread==8960 & zsh>0.590735
replace zpuma=3507 if metaread==8960 & zsh>0.6705906
replace zpuma=3508 if metaread==8960 & zsh>0.7685986
replace zpuma=3509 if metaread==8960 & zsh>0.8611774
replace zpuma=1300 if metaread==9040 & zsh>0
replace zpuma=1401 if metaread==9040 & zsh>0.3663274
replace zpuma=1402 if metaread==9040 & zsh>0.5648534
replace zpuma=1403 if metaread==9040 & zsh>0.7741507
replace zpuma=101 if metaread==9160 & zsh>0
replace zpuma=102 if metaread==9160 & zsh>0.284519
replace zpuma=103 if metaread==9160 & zsh>0.5445747
replace zpuma=104 if metaread==9160 & zsh>0.8010474
replace zpuma=2200 if metaread==9240 & zsh>0
replace zpuma=2300 if metaread==9240 & zsh>0.4063521
replace zpuma=901 if metaread==9260 & zsh>0
replace zpuma=902 if metaread==9260 & zsh>0.4260218
replace zpuma=3201 if metaread==9280 & zsh>0
replace zpuma=3202 if metaread==9280 & zsh>0.2833986
replace zpuma=3203 if metaread==9280 & zsh>0.7212811
replace zpuma=901 if metaread==9320 & zsh>0
replace zpuma=902 if metaread==9320 & zsh>0.2403697
replace zpuma=903 if metaread==9320 & zsh>0.4563083
replace zpuma=904 if metaread==9320 & zsh>0.6205314
replace zpuma=2300 if metaread==9320 & zsh>0.8098816
ren zpuma puma
collapse (sum) ct, by(ethn puma metaread) fast
egen cell=group(puma metaread)
egen msatot=sum(ct), by(metaread)
egen celltot=sum(ct), by(cell)
egen Emsatot=sum(ct), by(ethn metaread)
egen Ecelltot=sum(ct), by(ethn cell)
gen Esh=Ecelltot/Emsatot
gen Tsh=celltot/msatot
gen Edev=(Esh-Tsh)^2
collapse (sum) Edev ct, by(ethn metaread) fast
collapse (mean) Edev [fw=ct], by(ethn) fast
gen n=`i'
append using ./temp/residential-isolation-sim1
save ./temp/residential-isolation-sim1, replace

}

compress
save ./temp/residential-isolation-sim1, clear
su if n==0
su

log close